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ABSTRACT 

We present up to two decades of spectroscopic monitoring of the double-peaked broad Ha emission 
lines of seven broad-line radio galaxies. These extremely broad, double-peaked lines are successfully 
modeled by emission from gas in the outer accretion disk, and their profiles vary on timescales of 
months to years. This variability is a valuable probe for understanding the structure and dynamics of 
the accretion disk. We characterize the long-term broad-linc profile variability in detail, and compare 
the evolution of the line profiles with predictions from a few simple physical models for dynamical 
processes in the accretion disk. We find no evidence for persistent, periodic variability that would be 
indicative of a precessing elliptical disk, or a circular disk with a long-lived, single-armed spiral or 
warp. We do find transient, periodic variability on the orbital timescale, possibly resulting from shocks 
induced by tidal perturbations, and variability in the wings of the profile which suggests changes in 
the emissivity of the inner accretion flow. Dramatic but localized profile changes are observed during 
flares and high-states of broad-line flux. In 3C 332, these changes can be explained by a slow, smooth, 
secular change in disk illumination. In Arp I02B and 3C 390.3, however, a simple disk model cannot 
explain the asymmetric profile changes, suggesting a need for more sophisticated models. We also 
observe sharp, transient features that appear in the blue peak of the objects, which require a well- 
organized velocity field. 

Subject headings: accretion, accretion disks - galaxies: active - line: profiles 



1. INTRODUCTION 

In the standard model for active galactic nuclei 
(AGNs), the central engine is powered by an accretion 
disk feeding a supermassive black hole. The first direct 
dynamical evidence for the presence of the nuclear ac- 
cretion disk assumed to be powering AGNs emerged in a 
class of broad emission-line AGNs, known as "double- 
peaked emitters," with extremely broad (FWHM = 
10,000-20,000 km s~^) double-peaked Balmer emission 
lines (Eracleous & Halpern 1994) with line profiles rem- 
iniscent of the double-peaked lines known to originate 
in the accretion disks of cataclysmic variable stars (e.g., 
Horne & Marsh 1986). The detailed shape of the "disk- 
like" line profiles are well fitted by emission from pho- 
toionized gas in a circular or elliptical accretion disk 
around a central black hole (Chen et al. 1989; Chen & 
Halpern 1989; Eracleous et al. 1995). Ultraviolet spec- 
troscopy of several double-peaked emitters obtained with 
HST (Halpern et al. 1996; Eracleous et al. 2004) revealed 
that the Mg II A2800 line also has a double-peaked pro- 
file similar to the optical Ha and H/3 line profiles and 
appears to share the same origin as the Balmer lines, 
namely, a photoionized accretion disk. In 1/3 of the 
objects, the Lya line is weak and narrow, which is a 
prediction of the photoionized accretion disk model (Du- 
mont & CoUin-Souffrin 1990a,b,c; Rokaki et al. 1992). 
In the disk model, Lya photons in the dense inner re- 



gions of the accretion disk are trapped by scattering 
and destroyed by coUisional de-excitation, resulting in 
a suppressed emissivity of Lya in the inner disk where 
the rotational velocity of the gas is highest. In 2/3 of 
the objects, the higher-ionization UV lines and Lya are 
strong, with broad, single-peaked profiles, which is a pre- 
diction of the accretion disk -I- wind model, in which a 
highly ionized radiatively accelerated accretion disk wind 
produces broad, single-peaked profiles with logarithmic 
wings (Murray & Chiang 1997). 

Alternative models have been proposed for the origin 
of the double-peaked line emission, including: a bipo- 
lar outflow, binary black holes with separate broad line 
regions, an anisotropically illuminated spherical distri- 
bution of clouds, and bloated stars with cometary tails. 
These models have been progressively ruled out by incon- 
sistencies with observational results and by basic physi- 
cal arguments. Eracleous et al. (1997) were able to re- 
ject the binary supermassive black holes scenario favored 
by Gaskell (1983, 1996) for three double-peaked emitters 
(Arp 102B, 3C 390.3, 3C 332), based on long-term mon- 
itoring of their broad double-peaked Balmer lines (see 
also Halpern & Filippenko 1988). For > 20 yr their line 
profiles showed no evidence for the systematic drift of the 
radial velocity of the peaks expected for a binary system. 
The stationary peaks also rule out the possibility that 
the broad peaks are associated with individual clouds or 
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bloated stars, since the peaks would then move on the 
dynamical timescale of the broad line region, which is 
not observed. A collection of clouds on random orbits 
(Goad & Wanders 1996) is not physically plausible since 
the coUisional timescale between clouds is on the order 
of a few dynamical times of the broad line region, not 
allowing the configuration of clouds to last long enough 
to produce the persistent double-peaked emission. In 
at least one double-peaked emitter, 3C 390.3, a bipolar 
outflow has been definitively ruled out by reverberation 
mapping, which showed no time lag between the response 
of the wings of the light-echo of the broad Ha, H/3 line 
profiles, as well as the broad Lya and UV C IV line 
profiles, which would be indicative of a radial outflow 
(Dietrich et al. 1998; O'Brien et al. 1998; Shapavalova et 
al. 2001; Sergeev et al. 2002). 

1.1. Demography of Double- Peaked Emitters 

The first double-peaked emission line was discovered 
in the H/3 profile of the broad-line radio galaxy (BLRG) 
3C 390.3 (Lynds 1968; Burbidge & Burbidge' 1971; Os- 
terbrock, Koski, & Phillips 1975), and its large width 
and distinct red and blue peaks were attributed to mo- 
tion in an accretion disk (Oke 1987; Perez et al. 1988). 
Subsequent discoveries of double-peaked Ha lines in the 
BLRGs Arp 102B and 3C 332 had profiles that were suc- 
cessfully modeled as emission from photoionized gas in 
the outer regions of a circular relativistic Keplerian ac- 
cretion disk around a central black hole (Chen et al. 1989, 
Chen & Halpern 1989; Halpern 1990). Motivated by the 
fact that at that time the three known double-peaked 
emitters were hosted by BLRGs, Eracleous & Halpern 
(1994, 2003) conducted a spectroscopic survey of moder- 
ate redshift {z < 0.4) radio-loud AGN, and discovered 21 
new double-peaked broad emission lines (20% of the sam- 
ple), of which 13 of the profiles could be fitted with a sim- 
ple circular disk model. Careful starlight subtraction of 
ground-based spectra, and narrow-slit HST spectroscopy 
have revealed the presence of double-peaked emission 
lines in a large fraction of nearby low-ionization nu- 
clear emission line regions (LINERs), the most common 
type of AGN, including NGC 1097, (Storchi-Bergmann 
et al. 1993), M 81 (Bower et al. 1996), NGC 4203 (Shields 
et al. 2000), NGC 4450 (Ho et al. 2000), and NGC 
4579 (Barth et al. 2001). The Sloan Digital Sky Sur- 
vey (SDSS; York 2000) spectroscopic AGN sample (z < 
0.332) detected 116 additional double-peaked emitters 
which made up 4% of the entire sample. Of these, 76% 
are radio-quiet, 12% are classified as LINERs, 40% fall 
in the parameter space of a circular accretion disk pro- 
file, and 60% require some form of asymmetry in the disk 
(Strateva et al. 2003). 

1.2. Accretion Disk Model 

The circular relativistic Keplerian accretion disk model 
from Chen & Halpern (1989) is described by the following 
parameters: i, the inclination of the disk; and ^o, the 
inner and outer radii of the line-emitting portion of the 
disk in units of gravitational radii {rg= GM/c^); q, the 
emissivity power-law index (e oc ^~'); and a, the local 
broadening in km . The local broadening parameter 
is added to smooth out the peaks of the profile, which 
otherwise would be cuspy due to the sharp inner and 



outer boundary to the line-emitting region (see Chen et 
al. 1989). Another version of the disk model uses a con- 
tinuous emissivity power law with a break radius, instead 
of local broadening, to match the smooth peaks of the 
profile. In this model, there is no outer boundary to the 
line-emitting region, but simply a characteristic break 
radius at which the slope of the power-law changes. The 
emissivity power-law index depends on the geometry of 
the continuum source. At radii much larger than the size 
of the continuum source, the flux illuminating the disk 
drops off as r~^, and thus the emissivity, which is propor- 
tional to the illuminating flux (CoUin-Souffrin & Dumont 
1989), also has a power-law index, q, equal to 3 (Chen & 
Halpern 1989). In the disk model, double-peaked emis- 
sion lines will arise in the extreme case in which the line- 
emitting portion of the disk, £_i < ^ < ^o, is relatively 
small (i.e., a ring), resiflting in displaced broad peaks 
with a signature blueshifted peak higher than the red- 
shifted peak due to Doppler boosting, and a gravitational 
redshift of the wings of the line. The disk model produces 
an excellent fit to the peaks and wings of the broad Ha 
and H/3 lines, but there is often a central broad emission- 
line component at lower velocities that is not accounted 
for by the disk model. This excess emission may origi- 
nate from the "standard" BLR associated with the more 
common single-peaked broad emission lines. Eracleous 
et al. (1995) adapted the circular accretion disk model 
to elliptical disks in order to fit the profiles of double- 
peaked emitters with a red peak stronger than the blue 
peak (not possible in a circular disk, in which the blue 
peak is always Doppler boosted to be stronger than the 
red peak). This model introduces eccentricity and phase 
angle parameters to the circular model described above, 
and 5 then refers to the pericenter distance of the ellip- 
tical orbits. 

The symmetric, and often logarithmic, profiles of broad 

single-peaked optical and UV emission lines in standard 
broad-line AGNs can be reproduced with many velocity 
fields and geometric configurations of the broad-line re- 
gion (BLR), including both rotational and radial motion 
(see review by Mathews & Capriotti 1985). Because of 
this degeneracy, the kinematics and structure of the BLR 
remains enigmatic. Disks arc a promising candidate for 
the BLR structure (Shields 1977), and they are also in- 
voked to explain many other properties of AGNs. Rever- 
beration experiments (see review by Peterson 2001) mea- 
sure the light-crossing time of the BLR in nearby Seyfert 
1 galaxies as 1 2 It-wecks (c.g.,Maoz et al. 1990,1991; 
Netzer 1990), approximately the size expected for an ac- 
cretion disk around a supermassive black hole. A pho- 
toionized accretion disk can produce "standard" single- 
peaked broad emission lines if either the disk is large 
(Rokaki et al. 1992; Jackson, Penston, & Perez 1991), is 
oriented face-on (Corbin 1997), or emits a wind (Mur- 
ray & Chiang 1997). Given the frequency of "disk-like" 
emitters discovered in nearby galaxies, broad emission 
lines originating from the nuclear accretion disk may be 
a universal phenomenon in AGNs. 

1.3. Testing Models with Long-Term Profile Variability 

Spectroscopic monitoring of double-peaked emitters 

has revealed that a ubiquitous property of the double- 
peaked broad emission lines is variability of the profile 
shape on the timescale of months to years (Miller & Pe- 
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terson 1990; Romano et al. 1998). This slow, systematic 
variability of the line profile is on the timescale of dynam- 
ical changes in an accretion disk, and has been shown to 
be unrelated to the shorter timescale variability {^-^ days) 
seen in the overall flux in the line, due to reverberation of 
the variable ionizing continuum. Patterns in the variabil- 
ity of the broad Balmer lines are often a gradual change 
and reversal of the relative strengths of the blueshifted 
and redshifted peaks. Newman et al. (1997) successfully 
modeled the periodic oscillation in the relative heights 
of the blue and red peaks of the broad line in Arp 102B 
over 4 yr as an inhomogeneity, or bright spot, orbiting 
within the line-emitting accretion disk. From the orbital 
period of the bright spot (2.2 yr) and its velocity, they 
were able to solve for the mass of the black hole from 
the orbital radius of the spot (1.5 x 10^^ cm) and the 
inclination angle determined from the profile fit. This 
yielded a central black hole mass of 2.2;°;^ X 10^ Mq. 

Periodic variability of the red and blue peak strengths 
has also been attributed to a precessing elliptical disk 
(NGC 1097: Storchi-Bergmann et al. 1995, 1997) a pre- 
cessing single-armed spiral (3C 332, 3C 390.3: Gilbert 
et al. 1999; NGC 1097: Storchi-Bergmann et al. 2003), 
and a precessing warp in the disk. Elliptical disks and 
spiral waves have been detected in the accretion disks of 
cataclysmic variables (Steeghs, Harlaftis, & Horne 1997; 
Baptista & Catalan 2000), and a radiation induced warp 
has been detected in the large-scale disk of the AGN 
NGC 4258 (Maloney, Begelman, & Pringle 1996). Asym- 
metries in an AGN accretion disk, such as eccentricity, a 
spiral arm, or a warp, will be induced by the tidal effects 
of a binary companion in orbit around the primary black 
hole, or by the close passage of a massive star cluster 
(Chakrabarti & Wiita 1993, 1994). If the disk is self- 
gravitating, a spiral wave can emerge from a growing 
density perturbation, as is observed in the massive disks 
of young stellar objects (Adams, Ruden, & Shu 1989; 
Laughlin & Korchagin 1996). Elliptical disks have also 
been proposed to be a natural consequence of the tidal 
disruption of a star by a supermassive black hole, when 
the tidal disruption debris settles into bound eccentric or- 
bits around the central black hole (e.g., Rees 1988; Syer 
& Clarke 1992). The tidal disruption event scenario is 
a favored interpretation of the dramatic, sudden appear- 
ance of double-peaked emission lines in several LINERs 
(NGC 1097, M 81) and the BLRG Pictor A (Halpern & 
Eracleous 1994; Sulentic et al. 1995). Storchi-Bergmann 
et al. (2003) successfully fitted the detailed variability 
of NGC 1097 with a logarithmic one-armed spiral with 
a pattern period physically consistent with the central 
black hole mass estimated from the Mbh-c relation. Spi- 
ral waves are a physically desirable model. They can play 
an important role in accretion disks, because they pro- 
vide a mechanism for transporting angular momentum 
outward in the disk, allowing the gas to flow inwards 
towards the central black hole. 

Long-term profile variability is thus a useful tool for ex- 
tracting information about the structure and dynamics 
of the accretion disk most likely producing the double- 
peaked emission lines. In this paper, we present up to two 
decade's worth of spectroscopic observations of a sam- 
ple of double-peaked emitters with the goal of testing 
the robustness of the accretion disk model, and models 
for physical processes in the accretion disk as well as al- 



ternative models outside the disk paradigm, by directly 
comparing the observed evolution of the line profiles with 
the model predictions. Confirming or ruling out any of 
these models brings significant progress in developing a 
model for the origin and motion of the broad-line gas in 
the "special case" of double-peaked emitters, as well as 
in the standard broad-line regions of normal AGNs. 

2. SPECTROSCOPIC MONITORING PROGRAM 

2.1. Observations 

A long-term observing program was initiated in 1991 
to systematically take moderate resolution (4-6 A) opti- 
cal spectra of 20 known double-peaked emitters, selected 
for having apparent magnitudes bright enough to be ob- 
served with 2 m class telescopes in a reasonable amount 
of exposure time. The majority of the galaxies in our 
monitoring program are BLRGs from the Eracleous & 
Halpern (1994, 2003) survey of moderate-redshift radio- 
loud AGNs. Included are objects with clearly double- 
peaked Balmer lines as well as objects whose Balmer 
lines are very asymmetric with a displaced broad peak 
and shoulders (e.g., 3C 227 and Mkn 668). The spec- 
tra were obtained with the Lick 3m, the KPNO 2.1m, 
the MDM 2.4m, the CTIO 1.5m telescope, and recently 
with the 9.2m Hobby-Eberly telescope. Here we present 
the monitoring of seven BLRGs: Arp 102B, 3C 390.3, 
3C 332, PKS 0235-1-023, Mkn 668, 3C 227, and 3C 382. 
The remaining seven BLRGs are presented in Lewis & 
Eracleous (2006). The broad Ha profiles of 12 of the 14 
BLRGs in our monitoring program have been successfully 
fit at one time with a circular or elliptical disk profile in 
the literature. Here we present an elliptical disk profile 
fit for Mkn 668, leaving only 3C 227 with a profile that is 
too complex to fit with a single disk component. Table 
[1] lists the BLRGs presented in this paper, along with 
their apparent visual magnitude, redshift, and Galactic 
reddening. 

In order to sufficiently sample the long-term variabil- 
ity of the line profile, where significant changes typically 
occur on the timescale of months, when possible, 3-4 ob- 
servations of the objects were taken per year. The log 
of spectroscopic observations for the objects is listed in 
Tables [2H8l including the telescope and exposure time. 
Figure[l]shows a histogram of the temporal spectroscopic 
coverage of the BLRGs presented in this paper. 

2.2. Data Reduction 

The spectra obtained before 1999 were reduced with 
standard IRAF routines for bias correction, flat field- 
ing, sky subtraction, cosmic-ray rejection, and spectrum 
extraction. The calibration of one-dimensional spectra 
(wavelength calibration, atmospheric extinction and ab- 
sorption band correction, and flux calibration) was car- 
ried out with custom-made programs developed by our 
group (described in Eracleous & Halpern 1994). Spec- 
tra obtained after 1999 were reduced with IDL routines 
written for our observing program, which include scaling 
and shifting of the atmospheric bands for optimal cor- 
rection, and the removal of bad pixels based on outliers 
from the mean local spatial profile. Spectrophotometric 
standard stars (BD +17°4708, BD +26°2606, HD 19445, 
HD 84937: Oke & Gunn 1983; Fukugita et al. 1996) were 
observed to correct for atmospheric absorption features 
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and for construction of the sensitivity function. The ex- 
traction window was chosen to include 90% or 95% of 
the total enclosed flux, depending on the dominance of 
galaxy starlight. For our study we isolated the broad Ha 
line, since it is the strongest double-peaked emission line, 
and its high signal-to-noise makes it possible to study 
complex profile variability in detail. 

The continuum near the broad Ha line was modeled by 
the least-squares fit of the linear combination of a stellar 
template and a straight line (Eracleous & Halpern 1994). 
The stellar template was constructed from a scaled spec- 
trum of a normal galaxy [NGC 4339 (EO), NGC 5576 
(E3), NGC 7332 (SO pec sp), NGC 7457 (SA(rs)O), or 
UGC 555 (SO/a)], which was chosen from the fit that 
gave the smoothest residuals at the strong Na I D stel- 
lar absorption feature near 5900 A. Correct scaling of 
the stellar template is important for removal of the con- 
tribution of stellar absorption lines to the broad lines. 
However, the removal of Ha absorption may not be per- 
fect, since Na I D and Ha could come from different stars, 
and there might be some gas-phase absorption in Na I D. 
The spectra were corrected for Galactic extinction, and 
shifted to their rest wavelengths using the mean redshift 
measured from parabolic fits to the peaks of the strong 
narrow emission lines in the spectrum. 

In order to isolate the broad component of the Ha 
emission lines, the narrow emission lines superimposed 
on the broad Ha hne (narrow Ha, [O I] A6300, A6363, 
[N H] AA6548,6583, and [S H] AA6716,6731) were fitted 
with Voigt profiles plus a quadratic function approximat- 
ing the broad "continuum" below it and then subtracted 
from the spectrum. A Voigt profile fits the wings of the 
line profile more accurately than a single Gaussian. Due 
to the uncertainty in the shape of the broad line below 
the narrow lines, the following constraints were put on 
the Voigt profile fits: 1) the widths of the lines were fixed 
to the width of the isolated narrow [O I] A6300 line (or 
narrow Ha line for objects where [O I] is too weak), and 
2) the [O I] and [N H] forbidden-line doublets were fixed 
to have a total intensity ratio of 3:1. In the case of Mkn 
668, the multiple- Voigt profile fit to the Ha -|- [N II] 
narrow-line complex underestimates the flux in the nar- 
row Ha line for some of the epochs. In order to improve 
the fits, an additional constraint for [N II]A6583/Ha= 
0.97 was added, in order to match the narrow-line ratios 
of the best-fit profiles. The continuum-subtracted spec- 
tra were then normalized by the flux in the narrow [O I] 
A6300 (or narrow Ha) line for relative flux calibration. 
The narrow emission lines are used for relative flux c;al- 
ibration since the light-crossing time of the narrow-line 
region is long(~ 10^-10^ yr), and thus the narrow-line 
fluxes should not vary on the timescale of the monitor- 
ing program. Absolute flux calibration is not critical to 
this study, since we are most concerned with measuring 
the changes in the shape of the line profile that reflect 
dynamical changes in the broad-line gas and not the 
absolute flux of the broad line in response to the con- 
tinuum. Since the narrow-line region is predominantly 
unresolved in the spectra, the line ratios should be con- 
stant for all the observations, regardless of the slit-width 
or seeing conditions of each observation. The scatter in 
the narrow-line ratios therefore reflects the systematic 
error introduced by the narrow-line fitting, and is typ- 
ically on the order of 10%. The velocity scale for the 



spectra shifted to their rest wavelengths was determined 
from the relativistic Doppler formula: 
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The narrow-line-subtracted, narrow- line- flux- 
calibrated broad Ha profiles of the objects for all 
the spectra in our monitoring program are shown in the 
Appendix. In the case of 3C 382, the [O I] A6300 narrow 
line was too weak to use as a template for line fitting and 
flux calibration, and the scatter in the narrow Ha-f [N 
II] multiple- Voigt profile fit was too large to use the 
narrow Ha line as a reliable flux calibrator. Therefore, 
the proflles were neither narrow-line-subtracted nor 
flux-calibrated, and the profiles are shown normalized 
by the flux in the regions of the broad Ha line not 
contaminated by the narrow lines. 

2.3. Characterization of Variability 

We characterize in detail the variability of the line pro- 
files on timescalcs longer than the light crossing time of 
the broad- line region (t^t = tblr/c ^ weeks). Vari- 
ability on such long timescalcs (months to years) is un- 
likely to be the result of reverberation of a varialole ion- 
izing continuum (cf. Blandford & McKee 1982; Peterson 
1993). Rather, it is caused by changes in the structure, 
velocity field, or emissivity of the line-emitting gas. Si- 
multaneously to changes in the line profile there may also 
be slow changes in the integrated fiux of the broad lines 
caused by long-term variations in the flux of the ionizing 
continuum. 

Mean and root mean square (rms) profiles were cre- 
ated from the broad Ha profiles normalized by their to- 
tal broad-line fitix. The profiles are normalized by the 
broad-line flux in order to remove the first-order multi- 
plicative response of the profile shape due to changes in 
the line fiux. The rms profile thus measures deviations 
of the shape of the profile from the mean profile shape. 
Sharp features occur in the rms spectra at the location 
of the narrow emission lines due to systematic errors in 
the narrow-line removal. The velocities of the red and 
blue peaks of the mean spectrum were measured from 
the flux-weighted centroid of the top 10% of the peaks. 
The red and blue peaks of the individual proflles were 
characterized by measuring the first three moments of 
the top 10% of the peak: Mo, the average flux [P(R), 
P(B)]; Ml, the flux- weighted velocity centroid [V(R), 
V^(B)]; and when the signal-to-noise was sufficient, M2, 
the dispersion [fT(R), cr(B)]. The moments of the peaks 
were measured instead of using a Gaussian fit, in order 
to make no assumptions about the shape of the peaks, 
which are often sharp and asymmetric. Only the top 10% 
of the peaks were used to measure the moments of the 
peaks in order to avoid confusion with the neighboring 
broad profile shape. The wings of the profile were charac- 
terized by the full width at quarter-maximum (FWQM) 
and the shift at quarter-maximum (QMS) measured from 
the central velocity at quartca-maximum. The flux- 
weighted centroid at quarter- maximum (QMC), the total 
integrated flux, and the integrated flux in fixed velocity 
intervals were also measured. The statistical errors in 
the measurement of the average fiux, centroid, and ve- 
locity dispersion of the peaks were calculated from the 
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signal-to-noise in the continuum at a rest wavelength of 
6900-6930 A (this wavelength range was chosen so that 
it was covered by all the spectra in the study). The 
noise per pixel of the continuum in all the spectra is bet- 
ter than Upix ^ 0.01, which corresponds to a statistical 
error in the average flux and centroid of the peaks of 
less than 5%. Peak velocities and velocity dispersions 
are only given for peaks that do not have narrow lines 
superimposed on them, as the narrow-line subtraction 
is not accurate enough to reveal the true shape of the 
peak below it, and thus does not allow for a useful mea- 
sure of velocity and width. The 1 a error in measure- 
ments that include regions superimposed by the narrow 
lines, including the total broad-line flux, is calculated 
from adding in quadrature the la error in the flux of the 
least-squares Voigt profile fit and the systematic error in 
the relative flux calibration measured from the scatter in 
the narrow-line ratios (which dominates the error). In 
order to facilitate future studies of these objects, all of 
the measurements for the objects are reported in Tables 

E-M 

3. OBSERVED TRENDS 
3.1. Arp 102B 

The mean spectrum of Arp 102B (Figure [5^) has red- 
shifted and blueshifted peaks with V{K) — -1-5600 and 
V{Q) = —5200 km s~^, respectively. The rms spectrum 
(Figure [^b) indicates strong variability in the shape of 
the profile near the peaks of the mean spectrum, with 
the strongest variability in the wings of the profile. Fig- 
ure [3] shows the a) flux calibrated relative to narrow [O 
I]A 6300, b) FWQM, c) QMS, d) i^(R)/F(B), e) QMC, f) 
F(B), and g) cr(B) of the broad Ha line over the entire 
duration of the monitoring program. The F{K) / F(Q) 
ratio is measured from the average flux in the velocity 
ranges corresponding to the top 10% of the red peak 
(-1-4900 to -f6400 km s'^) and blue peak (-4500 to 
—5900 km s~^) of the mean spectrum. Fortunately, the 
velocity range for the red peak does not overlap with the 
location of the narrow [S II] lines. 

Flares in broad-line flux. - A rapid flare in the broad- 
line flux occurred in 1989, with a duration of ^ 10 days, 
and an increase in flux of 35%. In addition, there was a 
year-long period of a high state in flux in 1998, during 
which a large amplitude flare occurred with a duration of 
< 123 days, and an increase in flux of 60%. Both flares 
are also seen in the continuum flux near II/3, and thus 
are not the result of systematic effects in the relative 
flux calibration due to changes in seeing or slit width. 
The times of the flares and the high-state of flux are 
plotted with open circles in Figurc[3l During both flares, 
the FWQM decreased to < 1.65 x 10"^ km s'^ and the 
width of the blue peak narrowed to o'(B) = 400 km s~^ . 
During the rapid flare, the QMS and QMC shifted by ~ 
—500 km s""'^ , as a result of a sudden decrease in the 
strength of the red wing during the flare. 

Systematic variations in FWQM. - The FWQM varies 
in a step-like manner, with a rapid increase from 1.67 
X lO'' to 1.80 X 10" km s^^ from 1988.71 to 1989.50, 
and a drop to 1.75 x lO" km s^^ in 1993, and then a 
drop to 1.67 X 10** km s^^ in 1998. The changes in the 
FWQM reflect systematic changes in the wings of the 
profile. Along with changes in the width of the profile, 
the QMS also shifts from a mean of +1000 km s^^ to a 



mean of -1-600 km s^"'^ after 1997, revealing a weakening 
of the red wing relative to the blue wing. 

Red peak higher than blue peak. - The shape of the 
profile in 1983 is quite distinct from the later spectra, 
with a strong, broad red peak and i^(R)/i^(B) close to 
1. In fact, spectra of the broad Ha line reported by Miller 
& Peterson (1990) taken in 1982 have F(R)/F(B) > 1. 

Oscillation of F(R)/F(B). - F(R)/F(B) oscillates 
from 1991 to 1994, which was reported by Newman et 
al. (1997) (and measured from slightly different wave- 
length intervals), and this oscillation of i^(R)/i^(B) reap- 
pears in 1999. The least-squares fit of the F(R)/F(B) 
ratio to a sine curve with weights set to l/Cstat for the 
times of osciflation from 1990 - 1994 and 1999 - 2005 
results in a period of 2.073 ± 0.018 and 2.010 ± 0.006 
yr, respectively. Our period for the first oscillation is 
within the errors of Newman et al.'s (1997) measure- 
ment of 2.16 ± 0.07 yr. The oscillation of F(R)/F(B) 
is also evident in the QMC, which has the same pattern 
of oscillation. During the oscillation of F(R)/F(B), the 
blue peak has a minimum velocity during a maximum of 
F(R)/i^(B) in 1990.86 (when the blue peak is weak), and 
a maximum velocity during a minimum of -F(R)/-F(B) 
in 1991.47 (when the blue peak is strong). The width of 
the blue peak spikes during the i^(R)/F(B) oscillation at 
times when the very top of the blue peak appears to be 
suppressed in strength, and does not always occur during 
a maximum of i^(R)/i^(B). 

Systematic drift of the blue peak velocity. - The veloc- 
ity of the blue peak varies by a maximum amplitude of 
500 km s^^ from the mean. The velocity systematically 
decreases from 1994.87 to 1996.31, from -5300 to -4800 
km s"""^, and then increases back to —5300 km s~^ in 
1996.69. 

Systematic narrowing of the blue peak. - The width of 
the blue peak systematically narrows to a minimum of 
400 km s~i from 1994.46 to 1994.75, when the shape of 
the blue peak is visibly sharper. 

3.2. 3C 390.3 

The mean spectrum of 3C 390.3 (Figure [4^) has a 
strong blue peak with V{Yi) = —3300 km s~^ and a 
red peak on a fiat shelf, with ^(R) = -1-4300 km s~^ . 
The rms profile (Figure [4)d) shows five distinct regions 
of strong profile variability in the red central part of the 
profile, peaking at a velocity of +1300 km s^^ , near 
the velocities of the red peak and the blue peak of the 
mean profile, and in the wings of the profile. It should be 
noted that there is a spike in the rms profile at —2500 km 
s~^ that corresponds to the strong _B-band atmospheric 
absorption feature at 6870 A, which is not always di- 
vided out perfectly by the atmospheric standard stars 
in the reduction of the spectra. Figure [5] shows the a) 
flux cahbrated relative to narrow [O I] A6300, b) FWQM, 
c) QMS, d) F(R)/F(B), e) QMC, and f) V{K)-V{B) of 
the broad Ha line over the entire duration of the moni- 
toring program. Figure [6] shows the a) V'(B), b) V'(R), 
c) cr(R), and d) cr(B) of the peaks of the profiles. The 
F(R)/F(B) ratio is measured from the average flux in 
the velocity ranges corresponding to the top 10% of the 
red peak (+3300 to +5900 km s'^) and blue peak (-1700 
to —4500 km s^^) of the mean spectrum. Figure [7] shows 
the integrated flux in the five coherent variable regions of 
the profile determined from the rms spectrum, relative to 
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the total flux of the broad hne: the blue wing (—4500 to 
-12,300 km s-i), blue peak (-1700 to -4500 km s^i), 
central peak (-1700 to +3300 km s"i), red peak (+3300 
to +5900 km s"^), and red wing (+5900 to +15,000 km 
s""'^). The dramatic changes in the shape of the profile 
are also demonstrated in the individual profiles plotted 
in comparison with the mean profile scaled to the total 
flux of the line (see Appendix). 

Figure E shows F{B)/F{R) defined by Veilleux & 
Zheng (1991) (hereafter VZ91) and the velocity of the 
blue peak over 40 yr, with data points from 1974 to 1990 
from VZ91 shown with open circles, and data from our 
monitoring program shown with filled circles. The defi- 
nition of F(R) from VZ91 is actually centered on the red 
centra/ peak located at +1800 km s^^ , and F{B) includes 
the blue peak and blue wing of the profile. The quasi- 
periodic oscillation that they reported does not continue 
after 1990, and the F(B)/F(R) ratio continues to rise 
due to the decreasing strength of the central peak (see 
Figure [7|). The linear decrease of the velocity of the blue 
peak from -4700 km 8"^ in 1968 to -3200 km s'^ in 
1988 was reported by Gaskell (1996), and used to fit a 
double-line spectroscopic binary model with a period of 
300 yr and a binary mass of 6.6 x 10^ M©. The linear 
trend in the blue peak velocity does not continue after 
1988, and the velocity levels off at —3200 km s^-^ , which 
was also noted by Shapovalova et al. (2001) from their 
monitoring of the blue peak velocity of the H/3 profile 
from 1995 to 2000. 

Shapovalova et al. (2001) reported an anticorrelation 
between the peak velocity separation in the H/3 pro- 
file (measured from difference profiles, in which the 
minimum-activity average spectrum was subtracted from 
spectra not normalized by their broad-line flux) and the 
total broad line flux, as well as an anticorrelation with 
the flux measured in the continuum at 5125 A. We do not 
see such a relationship between peak velocity separation 
and total flux in our monitoring of Ha; however we mea- 
sure the peak velocities from the narrow-line-subtracted 
spectra, and have peak velocity separations that range 
from (6-8) x 10^ km s~^, while their measurements of 
V^(R)-V^(B) are in the range (8 - 15) x 10^ km s^^. It 
seems that their measurement of the peak velocity sepa- 
ration is affected when producing the difference profiles. 

Fluctuation of F (R)/F (B). - A fiuctuation of 
F{R)/F{B) is observed from 1994 to 1998 over 4 yr, 
with an amplitude of 0.11. During the maximum of 
F{K)/F{B) in 1996.12, the velocity of the blue peak 
drops by +600 km s~^ , and the width of the blue peak 
jumps by 90%. After the maximum of F(R)/F(B), the 
width of the blue peak steadily decreases, and the blue 
peak velocity increases. 

Flattening of the red peak. - The dramatic changes in 
the width and velocity of the red peak starting in 1997 
are a consequence of the shape of the red peak flattening 
and the red wing forming into a sharp shelf in 2000 (see 
individual spectra in the Appendix). 

Blueward drift of QMS and QMC. - After reaching a 
peak velocity in 1995, the QMS and QMC drift by ~ 
— 1000 km s~^ , reflecting the increase in strength of the 
blue peak and blue wing of the profile relative to the red 
wing. 

Transient sharp features in the blue peak. ~ Transient 
sharp spikes in the blue peak (shown in Figure [9|) appear 



in 1989.50, persisting for > 2 yr with a FWHM ~ 750 km 
s^^ and a velocity of -4300 km s~^, and again in 1998.27 
with a FWHM ~ 700 km s'^ that drifts from -3000 to 
—2600 km s^^ over a period of 1.5 months. Although the 
velocity of the spike in 1998 is quite close to the location 
of the atmospheric _B-band feature at —2500 km s^^, it 
is a real feature that is also detected in Shapovalova et 
al. 's (2001) monitoring of the H/3 profile of 3C 390.3 
taken with a different telescope from 1998.15 to 1998.54. 

Stronger blue peak during high state of flux. - The 
total flux of the broad line is in a high state in 2004 and 
the profile has a stronger and sharper blue peak than in 
the mean profile. However, it also appears that the blue 
peak is getting systematically sharper and stronger after 
2000 (Figure |42| . thus this behavior may be unrelated to 
the increase in flux. 

3.3. 3C 332 

The mean spectrum of 3C 332 (Figure [TUk ) has well- 
separated red and blue peaks with V(R) — +8300 km 
s~^ and V{B) — —6400 km s~^, respectively. The rms 
spectrum (Figure fTOb) is strongest in the red wing, and 
also shows profile variability in the blue wing. Figure [TT] 
shows the a) fiux calibrated relative to narrow Ha, b) 
FWHM, c) HMS, d) F(R)/F(B), e) HMC, and f) V^(R)- 
V^(B) of the broad Ha line over the entire duration of the 
monitoring program. Figure fT2l shows the velocity of the 
red and blue peaks of the profiles. The width, shift, and 
centroid is measured at half-maximum for this object, 
since the quarter-maximum is often located too far out 
in the wings of the spectrum to be measured reliably. 
The F(R)/F(B) ratio is measured from the average flux 
in the velocity ranges corresponding to the top 10% of 
the red peak (+7200 to +10,000 km s^^) and blue peak 
(—5300 to —8000 km s~^) of the mean spectrum. 

Maximum in F (R)/F (B). - F(K) /F(B) rises to a max- 
imum in 1995.43, at which point the velocity of the blue 
peak shifts to a lower velocity, and the velocity of the red 
peak shifts to a higher velocity. During this maximum, 
the HMC rises to a maximum red velocity of +2500 km 
s~^ , and then returns back to +500 km s"-'^ in 1998, 
where it remains constant. 

Changes in profile shape during high state of flux. - 
The FWHM and V{R)-V{B) begin to decrease when the 
flux increases in 1998. Figure [13] shows strong negative 
correlations between flux and the peak velocity separa- 
tion, FWHM, and the integrated flux in the red wing 
(+10,000 to +12,800 km s'^ ) and the blue wing (-8000 
to —10,800 km s~^ ), normalized by the total flux. These 
negative correlations indicate that at higher fluxes, the 
peaks shift closer together, and the wings of the profile 
are weaker. There also appears to be a bluer HMS dur- 
ing the high-state of flux, indicating a weakening of the 
red wing relative to the blue wing. 

3.4. PKS 0235+023 

The mean spectrum of PKS 0235+023 (Figure [Hi) 
has a narrow blue peak with V^(B) = —4700 km s^^ , 
and red peak with ^(R) = +2700 km s"^ . The rms 
profile (Figure [Hb) shows variability in the shape of the 
profile blueward of the blue peak, and redward of the red 
peak of the mean spectrum. Figure [TSl shows the a) flux 
calibrated relative to narrow Ha, b) FWQM, c) QMS, d) 
i^(R)/i^(B), e) QMC, and f) V{B) of the broad Ha line 
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over the entire duration of the monitoring program. The 
F(R)/i^(B) ratio is defined from the integrated flux in 
the regions of maximum variability in the rms spectrum 
near the red peak (+2300 to +5200 km s^^ ) and blue 
peak (—3900 to —6900 km s^^ ) of the mean spectrum. 

Dramatic change in profile shape after 1991. - The 
drop in F(R)/F(B) by 50%, the drift in the velocity of 
the blue peak by +500 km s""'^ , and the drift in velocity 

of the QMS and QMC by 2000 km s^^ , reveal a 

dramatic change in the profile shape between the 1991 
profile and the profiles after 1996, when the entire profile 
appears to get skewed towards the blue. 

Linear increase ofV(B). - The blue peak velocity in- 
creases roughly linearly from —3600 km s~^ in 1991 to 
-5000 km s~i in 2005. 

Minimum in FWQM. - The minimum in the FWQM 
of the profile from 1998-1999 is the consequence of a 
stronger blue peak as well as an almost complete disap- 
pearance of the red peak and red wing of the profile. 

3.5. Mkn 668 

The mean profile of Mkn 668 (Figure [T6k) has a strong 
central red peak with V^(R) — +1700 km s~^ , and a 
broad base with a blue shelf. The rms spectrum (Figure 
116b ) shows strong variability redward and blueward of 
the central red peak of the mean profile. Figure [T71 shows 
the a) flux calibrated relative to narrow Ha, b) FWQM, 
c) QMS, d) F(R)/F{B), e) QMC, and f) 1/(R) of the 
broad Ha line over the entire duration of the monitoring 
program. The F(R)/F(B) ratio is measured from the 
integrated flux in the regions of strong variability mea- 
sured by the rms spectrum, redward of the central peak 
(+3200 to +5600 km s~^ ) and blueward of the central 
peak (—800 to —3200 km s^^ ) of the mean profile. 

Fluctuation of F(R)/F(B). - The F(R)/F(B) ratio 
fluctuates, with a maximum of 2.5 in 1990.42, decreasing 
to a minimum of 1.5 from 1996 to 2001, and then in- 
creasing to a maximum of 2.3 in 2002. These changes in 
i^(R)/i^(B) are reflective of the dramatic changes in the 
shape and peak velocity of the central peak of the profile. 
The velocity of the central peak shifts from a maximum of 
+2500 km s^^ during the maximum of i^(R)/i^(B), and 
to a minimum of < +1700 km during the minimum of 
i^(R)/i^(B), while the shape of the profile changes from 
a skewed to a more symmetric shape. 

Dramatic changes in V (R). - The historical velocity 
curve of the central red peak is shown in Figure [TS] with 
data points in 1974 from Osterbrock & Cohen (1979), in 
1982 from Gaskell (1983), and data points from Marziani 
et al. (1993) from 1985 to 1991 measured from the H/3 
profile. The central peak was at a maximum velocity 
of +2700 km s~i from 1974 tol982, before shifting to a 
lower velocity of +1800 km s^^ in 1985. When the peak 
shifts to smaller velocities, the peak velocity measure- 
ment of the Ha profile is affected by systematic errors in 
the narrow-line subtraction. Due to the uncertainty in 
the exact shape of the broad-line underneath the narrow 
Ha+[N II] narrow- line complex, upper limits are shown 
for our measurement of V^(R) when the measured velocity 
is within the velocities of the narrow lines. 

3.6. 3C 227 

The mean profile of 3C 227 (Figure [T^ ) has a strong 
central blue peak with T^(B) = —1500 km s~^ , and a 



broad base with a red and blue shelf. The rms profile 
reveals a strong and narrow region of variability on the 
blue side of the central peak, with a FWHM ~ 2200 
km s^"'^ , and strong variability redward of the central 
peak, as well as some variability in the shape of the blue 
and red shelf of the profile. Figure [20] shows the a) flux 
calibrated relative to narrow Ha, b) FWQM, c) QMS, d) 
i^(R)/i^(B), e) QMC, and f) ^(B) of the broad Ha line 
over the entire duration of the monitoring program. The 
i^(R)/i^(B) ratio is measured from the integrated flux 
in the regions of strong variability measured by the rms 
spectrum redward of the central peak (+1700 to +3400 
km s^^ ) and blueward of the central peak (—1200 to 
—4400 km s^^ ) of the mean spectrum. 

Excess feature in the blue peak. - The Fiji)/ F{Q) ra- 
tio increases monotonically from 1990 to 2000, reveal- 
ing an excess feature in the central blue peak that fades 
away. The historical velocity curve for 3C 227 (Figure 
[2T|) . which includes data points in 1974 from Osterbrock, 
Koski & Phillips (1976), and in 1982 from GaskeU (1983), 
indicates that the excess in the central blue peak per- 
sisted from 1974 to 1998. The QMS and QMC increase 
by +800 km s"^ between 1995 and 2000, also a result of 
the fading away of the blue feature. 

3.7. 3C 382 

The mean profile of 3C 382 (Figurel^^) has a red peak 
and a blue shoulder that are superimposed by the Ha+[N 
II] narrow-lines. The mean and rms profiles are mea- 
sured without the narrow emission lines subtracted from 
the spectrum, because of the large errors in fitting the 
weak [O I] A 6300 narrow line and the crowded Ha+ [N II] 
narrow-line complex. The rms profile (Figure [22b ) shows 
variability in the shape of the red wing, as well as a nar- 
row region of profile variability at the velocity of the red 
peak of the mean profile, and in the far blue wing. Figure 
1231 shows the integrated flux in the regions of maximum 
variability measured by the rms spectrum, relative to 
the total integrated flux of the line, avoiding regions su- 
perimposed by narrow lines, near the red peak (+2100 to 
+3800 km s~i ), the red wing (+4400 to +6200, +8000 to 
+15,000 km s^i ), and the blue wing (-7500 to -11,300 
km s^^ ) of the mean spectrum, over the entire duration 
of the monitoring program. i^(R)/F(tot) drops by 15% 
in 1994 to 1995, and then returns back to its mean value 
in 1997. During the drop in i^(R)/i^(tot), the normalized 
flux in the red wing increases by 25%, and then steadily 
decreases. After this jump in F(RW)/F(tot), there is a 
sudden increase in F(BW)/F(tot) in 1998 by 40%, after 
which the normalized flux in the blue wing steadily de- 
creases. After 2002, while the normalized flux in the red 
and blue wing is decreasing, the normalized flux in the 
red peak increases. 

Fluctuation of F (RW)/F (BW).- The one cycle of fluc- 
tuation in the F(RW)/F(BW) ratio between 1994 and 
2000, is a result of the sharp increase in the relative 
strength of the red wing in 1995 followed by an in- 
crease in the blue wing in 1998. This fluctuation of 
i^(RW)/i^(BW) does not repeat during the 16 yr of our 
spectroscopic monitoring program, and thus if this trend 
is periodic, its timescale is > 20 yr. 

4. COMPARISON WITH MODELS 
4.1. Arp 102B 
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4.1.1. Circular Disk 

The characteristic timescale of the F(R)/i^(B) oscilla- 
tion in Arp 102B (~ 2 yr), is on a much shorter timescale 
than the slow, step-like variations in the FWQM (see Fig- 
ure [3]). In order to decouple the oscillation of the peaks 
of the profile from the slower, systematic changes in the 
wings, the step-like variation in the FWQM of the pro- 
files was successfully modeled by changing the inner and 
outer radii of the best-fit circular disk model from 1990- 
1994 {i = 31°, = 300, Co = 825, g = 3, cr = 1050 
km s^^). Figure [24] shows the three model disk profiles 
that fit the mean FWQMs of the profiles plotted with 
solid lines in Figure [3Jd. The step-like variations of the 
inner and outer radius (from C,; — 400 and — 700, to 
C, = 325 and = 825, to C» = 300 and Co = 825) of 
the disk are either a result of changes in the emissivity 
of the disk, perhaps due to changes in the illumination, 
or changes in the structure of the accretion flow itself. 
The radial structure of the accretion disk evolves on the 
viscous timescale, Tmsc ~ r'^ / {TdynCicl) , where Tdyn is 
the dynamical time, a (~ 0.1) is the Shakura -Sunyaev 
viscosity parameter (Shakura & Sunyaev 1973), and Cs is 
the sound speed. This timescale is several orders of mag- 
nitude larger than the time over which the apparent radii 
of the disk change (~ 10 yr); thus the structure of the 
disk cannot be varying fast enough to cause the changes 
in the radii of the line-emitting region of the disk. If il- 
lumination is the driver for the changes in radii, then it 
follows that since the broad-line flux does not correlate 
well with the changes in radii, it may be that the geom- 
etry, not the flux, of the continuum source is changing 
the emissivity profile of the disk. 

The profiles were fitted with one of the three model 
disk profiles with the closest FWQM, by scaling the 
model to minimize the sum of the squares of the errors 
from the model (least-squares fit). Regions of the profile 
contaminated by the narrow lines were not used in the 
fit. The scaled model disk profile for each spectrum is 
shown in blue in the spectra in the Appendix. During 
times of the oscillation of F(R)/F(B) the profiles were 
fitted with a model disk profile, allowing for a Gaussian 
excess above the disk profile. Our method is similar to 
that of Newman et al. (1997) for modeling the profiles 
with a disk plus a Gaussian, except that we keep our disk 
profile fixed at any given time, and allow the Gaussian 
to vary, while Newman et al. (1997) allowed both com- 
ponents to vary, by allowing the inner and outer radii 
of the disk model to change. Figure [^5] plots the peak 
velocity of the Gaussian excesses during the oscillations 
of F(R)/F(B) as a function of time. The velocity of the 
excesses are well constrained when they have a velocity 
near the peaks of the profile. However, there are large 
uncertainties when the excesses are located within the ve- 
locities of the narrow Ha-|-[N II] complex, since the shape 
of the broad line underneath the subtracted narrow lines 
is not well constrained due to errors in the narrow-line 
subtraction. The FWHM of the Gaussians fits from 1991 
to 1994 are ~ 4000 km s~^, while the Gaussians are nar- 
rower from 1999 to 2003 with a FWHM of ~ 2000 km 
s~^, and narrow even further to 1000 km s^^ in 2004. 
The most natural explanation for an excess that peri- 
odically oscillates between the red and blue sides of the 
profile is rotation. 



4.1.2. Orbiting Bright Spot? 

Newman et al. (1997) fitted the first episode of the 
i^(R)/i^(B) oscillation from 1991 to 1995 with a bright 
spot at a radius of 455 r^, with an angular width of 8° 
(FWHM = 2000 km s-i), a radial width of 5 r^ and an 
orbital period of 2.2 yr. Newman et al. (1997) determined 
the radius of the bright spot from the peak velocity of 
the Gaussian fit to the excess in the blue peak in the 
1991 June 17-20 spectra. If the bright spot is assumed 
to be at a phase angle (0) of 270° in the disk (when the 
bright spot is at its maximum blue radial velocity), the 
gravitational radius of the bright spot can be calculated 
from the relativistic Doppler factor, D = v/ve, which in 
the weak-field approximation is 



D={1- 3/0^/^(1 + C'''^sini sincj))-', (2) 

where i is the inclination determined from the disk fit, 
and the observed velocity is then 



-1/2, 



D 



1 



(3) 



The radius, period, and phase, = (27r/P)(t — ^0)7 of the 
bright spot for each time of oscillation of i^(R)/i^(B) is 
determined from a least-squares fit of the observed veloc- 
ity curve of the Gaussian excesses to the velocity curve 
of a bright spot in the disk (eq. [3]) , with weights set to 
1/cr^. The resulting radii for the first and second episode 
of the oscillation are 427 ± 5 and 757 ± 4 respec- 
tively, and periods of 2.247 ± 0.003 and 2.031 ± 0.001 
yr, with quoted errors equal to 1 a. The quoted errors are 
the standard deviation of the best-fit parameters for the 
model velocity curve and do not reflect systematic errors 
from the assumption that the data behave like the model. 
Combining the orbital periods of the bright spots with 
their radius in units of r^, a dynamical measurement of 
the central black hole can be made. Using the complete 
form of Kepler's third law for a circular orbit. 



P = 2ti 



GMbh 



1/2 2TrGMf/^ 



(4) 



this yields central black hole masses of 2.60 ± 0.05 xlO* 
M0and 1.00 ± 0.01 x 10^ Mq. Our measurement of the 
black hole mass from the first oscillating bright spot is 
consistent with Newman et al. 's (1997) measurement; 
however, it is discouraging that modeling the second 
episode of the F(R)/F(B) oscillation with an orbiting 
bright spot yields a discrepant black hole mass. An inde- 
pendent estimate of Arp 102B's central black hole mass 
was made by Lewis & Eracleous (2006) from the mea- 
sured stellar velocity dispersion (Barth et al. 2002), and 
the correlation between the velocity dispersion and the 
mass of the central black hole (Tremaine et al. 2002). 
Including the uncertainty in the slope of the Mbh-c 
relation, they determined a central black hole mass of 
(1.1 ± 0.2)xl0^ Mq. It may be that a single orbit- 
ing bright spot is not a correct interpretation of the 
Gaussian excesses, and the excesses may be the conse- 
quence of a more complex configuration of emissivity 
enhancements orbiting in the disk. The transient na- 
ture of the oscillations suggests a time-dependent con- 
figuration. One mechanism to produce transient bright 
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spots in the disk is spiral shocks induced by tidal per- 
turbations. Chakrabarti & Wiita (1993) conducted two- 
dimensional hydrodynamical simulations to investigate 
the tidal effects on an AGN accretion disk by passing 
massive perturbers ranging from Mpert/MB}i= 10^^ - 
10~^. Stellar-mass perturbers, with Mpert/^^BH< 10^^, 
were not included in the simulation due to limitations 
in their grid resolution. The simulations revealed den- 
sity waves in the disk that steepened into spiral arms 
and knotty structures of increased density, temperature, 
and emissivity, that formed, fragmented, and reformed. 
The lifetimes of these features was longer than the or- 
bital period in the inner regions of the disk, and during 
their lifetime, they propagated outward in the disk. This 
scenario may be a better physical model for the "bright 
spot" in Arp 102B, which appears for two orbital periods 
with a radius of 427 r^and then reappears several orbital 
periods later, further out in the disk at a radius of 757 

4.1.3. Collection of Orbiting Clouds in a Disk? 

An alternative model that predicts an oscillation of 
F{R)/F{B) is not a single orbiting bright spot, but 
rather a disk-like distribution of clouds in nearly Kep- 
lerian orbits, with randomly distributed clouds having 
enhanced emissivity (Sergeev et al. 2000, hereafter SOO). 
Monte Carlo simulations by SOO of clouds with a distribu- 
tion of emissivities produce times of sinusoidal variation 
of the red and blue peaks, as well as a sinusoidal variation 
of the central portion of the profile with half the period 
of the F(R)/F{B) variation (since an orbiting cloud will 
pass through the center twice per orbit). Contrary to 
this, SOO's measurement of the central-to-total flux ratio 
in their spectra of Ha from 1992 to 1995 had a sinusoidal 
variation with the same period within the errors to that 
of the F{R)/F{B) oscillation. In order to explain this 
discrepancy with the predictions of the model, they pos- 
tulated that the line emission originates in the face of 
the cloud illuminated by the central continuum source. 
Thus, if the disk is edge-on, clouds traversing in front of 
the continuum source will have their emitting face hidden 
from the line of sight, and their contribution to the line 
flux is seen to traverse the central portion of the profile 
only once per orbit. 

Figure [26] shows the red, blue, and central flux ratios 
defined by SOO, measured for the spectra in our moni- 
toring program. Plotted with solid lines are the best- 
fit sine curves for their data set from 1992 to 1995. 
These curves are extended with dotted lines over the 
baseline of our observing program. It should be noted 
that these fiux ratios include regions superimposed by 
narrow lines, and thus are subject to systematic errors 
introduced by the narrow-line subtraction, and typical 
error bars are shown. Our data fit their model for 
all the ratios from 1991 to 1994 [during the first os- 
cillation of F(R) / F(B)], and follow reasonably well for 
i^(R)/i^(tot), F(B)/F(tot), and F(R)/F(B) during the 
later oscfllation of F(R)/F(B) between 1999 and 2004. 
Contrary to their model, the F(C)/i^(tot) ratio appears 
to stop oscillating after 2001, although the error bars in 
the F{C) /F{tot) ratio are too large to make a convincing 
argument either way. 

SOO simulated the residuals from the mean normalized 
profile with a collection of 1500 clouds with an inter- 



nal velocity dispersion resulting in a FWHM of 2000 km 
s~^. Their observed residuals from the scaled mean pro- 
file show enhancements traversing from the red to the 
blue side of the profile only, with no tracks from the blue 
to the red side. They again attribute this to the phase ef- 
fects that explain the equal periods of i^(C)/F(tot) and 
F(R)/F(B). Due to the uncertainty in the narrow-line 
subtraction in the center of the profile, as well as the un- 
known nature and variability of the central broad com- 
ponent of the profile not accounted for by the disk-fit (see 
§??), it is not clear whether a reliable detection of tracks 
traversing through the central part of the profile is possi- 
ble. Our measurement of F(C)/i^(tot) and the residuals 
from the disk profile do not show definitive evidence for 
phase effects, and thus it seems that the edge-on con- 
figuration does not seem necessary, in either the cloud 
model or the bright spot model (in which the bright spot 
is also expected to traverse the center twice per orbit). 
In general the cloud model is problematic, because the 
clouds must have large internal velocity dispersions and 
a large covering factor, which would lead to dissipation of 
the clouds or destruction by collisions on the timescale of 
decades (SOO), unless there was a mechanism to produce 
new clouds on this timescale. Thus, without the neces- 
sity for phase effects, the orbiting cloud model seems less 
compelling than emissivity enhancements orbiting in an 
accretion disk. 

4.1.4. Precessing Spiral Arm, Eccentric Disk, or Warp? 

The bright spot model is not the only variation of the 
disk model that can produce an oscillation in F(R)/F(B) 
and an excess above the disk profile that shifts from one 
peak to the other periodically. A precessing eccentric 
disk, and a disk with a precessing spiral arm or warp, 
will have periodic changes of F(R)/F{B). In the eccen- 
tric disk model, the precession of the disk causes a si- 
nusoidal variation of F(R)/F(B) as well as a sinusoidal 
variation of the peak velocities with the same period, 
with a minimum blue peak velocity during the minimum 
of i^(R)/i^(B) and a minimum red peak velocity during 
the maximum of F(R)/F(B) (for emissivity power-law 
index q < 3; see Lewis et al. 2004). The lack of a shift 
in the velocity of the peaks on the same timescale as the 
i^(R)/i^(B) oscfllation in our data does not support this 
model. Most importantly, the timescale of the precession 
of the disk is determined by either the general relativis- 
tic advance of pericenter or the tidal effects of a binary 
companion, and is on the timescale of a few centuries for 
a 10* Mq black hole (Eracleous et al. 1995). The 2 yr 
period of the F(R)/F{B) oscillation is much too short 
to be associated with a precessing disk or a precessing 
warp (see Storchi-Bergmann et al. 1997, eq. [4]), and its 
transience could not be associated with a persistent phe- 
nomenon. 

The timescale of the F(R)/F(B) oscillation due to a 
precessing spiral arm is determined by the pattern pre- 
cession period, which is typically an order of magnitude 
longer than the dynamical time, and less than the sound- 
crossing time (see discussion in Storchi Bergmann et 
al. 2003, and references therein). For a 10* Mq black 
hole and an inner radius of 300 r^, that translates to 
a period between 1 and 20 yr (using eqs. [5] and [6] in 
Storchi-Bergmann et al. 2003). The 2 yr period of the 
i^(R)/i^(B) oscillation is within the range of plausible 
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periods for a spiral arm; however, another consequence 
of a spiral arm are changes in the velocity and width of 
the peaks of the profile as the arm precesses. For a one- 
armed spiral that spans the entire radial extent of the 
line-emitting disk, during the minima of F(R)/F(B) the 
blue peak narrows and is at a maximum velocity, while 
during the maxima of F{IC)/ F(B) the red peak narrows 
and is at a maximum velocity (see Gilbert et al. 1999, 
Storchi-Bergmann et al. 2003, and Lewis et al. 2004 for 
examples). In our data, the width of the blue peak does 
not appear to narrow during minima of the F(R)/ F{Ji) 
oscillation, nor does the velocity of the blue peak in- 
crease, except during the minimum of F(R)/F(B) in 
1991.4. Thus, with the only spiral-like behavior being 
the periodic variation of F(R)/F{B), there are no com- 
pelling reasons to favor a spiral arm over a bright spot or 
other configuration of emissivity enhancements orbiting 
in the disk. 

4.1.5. Reverberation of Flares in the Continuum? 

The change in shape of the line profile during the small 
amplitude, rapid flare in 1989 July (Figure 2.13), is on 
a short-enough timescale to be the result of light travel 
time delays of the reverberation of the variable ionizing 
continuum. Using the black hole mass estimated from 
two independent methods: the bright spot model and 
the stellar velocity dispersion (Lewis & Eracleous 2006) 
which both yield a mass of ~ 1 x 10* M©, and the outer 
radius of the best-fit disk profile in 1989 of 825 r^, the es- 
timated light-travel time across the broad-line region is 5 
days. A flare in the continuum with a duration on the or- 
der of the light-crossing time and much greater than the 
recombination time [rfiare ~ tlt >> Tree = ('T-etts)"^, 
where Ue is the electron density and as is the case B 
recombination coefficient], can cause a detectable light 
echo to propogate across the profile (Peterson 1988). A 
light echo from such a flare is expected to produce spikes 
in the profile that drift from higher to lower velocities 
as it propagates outward in the disk (Stella 1990). Dur- 
ing the flare in 1989, the line profile has a dramatically 
weaker red wing (Figure [27|) , and a narrower and stronger 
blue peak. The flare in 1998 occurs over a much longer 
timescale than the light-crossing time, and appears to 
have the same changes in the profile shape, namely a 
narrower red wing and stronger blue peak. Thus light- 
travel time effects are unlikely to be causing the dramatic 
changes in the profile shape during both flares. 

During a flare, the increase in central continuum emis- 
sion can temporarily illuminate the regions further out 
in the disk, increasing the contribution of emission from 
lower velocities. Increasing the outer radius of the circu- 
lar disk model mainly broadens the peaks of the profile 
and shifts them to lower velocities but does not signifi- 
cantly affect the wings of the profile. Alternatively, the 
rapid flare could temporarily fully ionize the inner regions 
of the disk, causing an increase in the effective inner ra- 
dius of the line-emitting portion of the disk or a flatter 
emissivity law at small radii. The result is a decrease in 
the strength of the wings, a shift of the peaks to lower ve- 
locities, and an increase in F(R)/F(B) (shown in Figure 
[24b). The decrease of F(R)/F(B) during the flares is not 
compatible with this scenario. In general, the decrease in 
F{R) /F{B) and in the FWQM cannot be modeled by in- 
creasing the central continuum illumination of a circular 



disk. 

4.1.6. Gravitational Lensing? 

Gravitational lensing can amplify a broad emission-line 
if the Einstein radius of the microlens is of a compara- 
ble size to the accretion disk producing the line emis- 
sion. An interesting consequence of this, is that the pro- 
file can be amplified asymmetrically if the lens is off- 
center. Popovic et al. (2001) calculated two scenarios in 
which a microlens could produce sinusoidal variations in 
F{K)/F{B) similar to those observed in Arp 102B: (1) 
a supermassive binary companion at a distance of 0.3- 
Ipc, which takes a few years to transit in front of the 
disk, but has an orbital period of ^ 1000 yr; (2) a 1 
AIq star in an intervening galaxy (l/20th the distance 
of Arp 102B) with a relative velocity of 200 km s~^ . 
The important difference between these scenarios is the 
frequency with which the lensing events would occur. In 
the supermassive binary case, the lens should only tran- 
sit in front of the disk once every 1000 yr, and thus only 
one period of the F(R)/F(B) sinusoidal variation should 
be observed during the entirety of our monitoring pro- 
gram. For the stellar-mass microlens case, the lensing 
transit events would randomly cause one period of an 
oscillation of i^(R)/i^(B) during the transit, with a fre- 
quency depending on the density of microlenses in the 
intervening galaxy. The multiple consecutive periods of 
the i^(R)/i^(B) oscillation observed after 1999 rule out a 
transiting scenario for the source of the oscillation. How- 
ever, microlensing could be a mechanism to amplify the 
line profile asymmetrically and cause a flare of broad- 
line flux. This could be a possible explanation of the 
stronger blue peak of the profiles during the flares, but 
the decrease in flux of the red wing during the flares is not 
consistent with this picture, nor is the coincidence of both 
flares having a stronger blue peak, which in the gravita- 
tional lens model, should be a random consequence of 
the alignment of the microlens with respect to the disk. 

4.1.7. Drift of Peaks due to Changing Outer Radius? 

During the systematic drift of the blue peak velocity 
from -5300 to -4800 km s^^ and back to -5300 km 
s~^ again from 1995 to 1996, it is difficult to determine 
whether the velocity of the red peak is also shifting, be- 
cause of the [S II] narrow lines superimposed on the red 
peak. In the circular disk model, the peak velocities will 
shift when the outer radius of the disk increases. If both 
of the peaks were drifting to lower velocities, then this 
could be modeled as an increase in the outer radius of 
the line emitting portion of the disk, which causes the 
peak separation to decrease. Figure [28] shows the spec- 
tra during the drift of the blue peak from 1995 to 1996 in 
comparison to the change in the profile shape caused by 
an increase in outer radius from 825 to 1300 r^. Although 
the shift in the blue peak velocity in the model matches 
the data, the shift of the red peak and broadening of 
both peaks in the model do not. This drift in the blue 
peak velocity occurs during a time of low i^(R)/i^(B). 
In general, there appears to be a trend towards a lower 
blue peak velocity when the blue peak is stronger, as is 
the case during flares in the broad-line flux. 

4.2. 3C 390.3 
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4.2.1. Binary Broad-Line Region? 

Eracleous et al. (1997) created yearly averaged velocity 
curves of the red and blue peaks of 3C 390.3 through 
1996 and found that the red peak velocities were not 
consistent with the motion of a spectroscopic binary, and 
the best-fit model yielded a lower limit on the binary 
mass of > 2.6 xlO^^ Mq, larger than any known central 
black hole mass, and over 2 orders of magnitude larger 
than the black hole mass inferred from the stellar velocity 
dispersion, Mbh= 5 ± 1 x 10^ Mq (Lewis & Eracleous 
2006). Eracleous et al. (1997) noted that the velocity of 
the red and blue peaks are highly sensitive to changes in 
the shape of the peaks, which is evident in our spectra of 
the transient sharp features in the blue peak in 1989 and 
1998 and in the flattening of the shape of the red peak 
from 1997 onwards, and occur on much shorter timescales 
than that predicted by an orbiting binary model. The 
reverberation mapping results for 3C 390.3 (see §??) also 
contradict the binary broad-line region scenario for the 
source of the displaced red and blue peaks of the profile, 
since no lag is detected in the response of the red and blue 
peaks, which would be expected if they were produced 
by a binary separated by 1.5 pc, as is determined from 
the best-fit binary model from Eracleous et al. (1997). 
It is also noteworthy that the widths of the two peaks 
of the Ha line are narrower than their separation, which 
is opposite from what one expects for a binary broad- 
line region governed by Kepler's laws (M. V. Penston 
1998, private communication; see footnote 3 in Chen et 
al. 1989). 

4.2.2. Bipolar Outflow or Orbiting Bright Spot? 

In general, it is difficult to fit the low F(R)/F(B), nar- 
row FWQM, and strong central component of the broad 
Ha profile of 3C 390.3 with a simple circular disk. The 
best-fitting disk model to the broad Ha line by Eracleous 
& Halpern (1994) in a spectrum obtained on 1988 July 
9 successfully fits the wings of the profile but underes- 
timates the strength of the blue peak. Even more diffi- 
cult, is reconciling the circular disk model with the 1976 
and 1980 Ha profiles reported in Zheng et al. (1991), in 
which F(R)/F(B) changed from 0.8 to 1.5, indicating a 
stronger red peak than the blue peak in 1980. Zheng et 
al. (1991) modeled the two epochs of the Ha profile with 
an edge-on (non-relativistic) circular disk plus a bright 
spot with an orbital period of 10.4 yr, and also with a 
double-stream outflow. In the bipolar outflow model, the 
changes in F(R)/F(B) are due to variations of the rela- 
tive intensity of the ionizing continuum illuminating the 
two cones, or due to changes in orientation of the cones, 
perhaps resulting from a precession of the rotation axis 
caused by its misalignment with the axis of the accretion 
disk. However, the theoretical timescale for the align- 
ment of the disk (and hence the outflow axis) is much 
too long, 10^ — 10^ yr, to account for the timescales 
observed (Natarajan & Pringle 1998; Natarajan & Ar- 
mitage 1999). Regular periodicity of i^(R)/i^(B) on the 
timescale of 10 yr is not evident in our monitoring data, 
ruling out the possibility of regularly precessing jets. The 
bipolar jet model is disfavored by reverberation mapping 
campaigns of 3C 390.3 starting in 1992, which show no 
lag between the response of the red side of the profile 
compared to the blue side in the Ha and H/3 lines, which 



would be observed in the case of radial motion (Dietrich 
et al. 1998; Shapovalova et al. 2001; Sergeev et al. 2002). 
Thus, a bipolar outflow, if present from 1976 to 1980, 
must have disappeared by 1992 and would not be able to 
explain the fluctuation of i^(R)/i^(B) from 1994 to 1998. 
The mass of the central black hole can be estimated by 
combining the light crossing time of the broad- line region 
of 20 days (Dietrich et al. 1998) with the best-fit outer 
radius of the accretion disk model fit of 1300 rg from Er- 
acleous & Halpern (1994), yielding a black hole mass of 
3 x 10^ Mq, which is in relatively good agreement with 
the black hole mass derived from the stellar velocity dis- 
persion. For a central black hole mass of 3 x 10* Mq, 
and an inner radius of 380 rg, this translates to orbital 
periods ranging from 2 to 14 yr. The two episodes of 
an oscillation of F(R)/F(B) in 1976 to 1980, and 1994 
to 1998, have periods that are both within the range of 
orbital periods in the disk and could be attributed to 
transient emissivity enhancements orbiting in the disk. 

4.2.3. Precessing Eccentric Disk, Spiral Arm, or Warp? 

The transient and irregular nature of the i^(R)/i^(B) 

oscillations disfavors a precessing eccentric disk, for 
which the period of the oscillation docs not change, and 
the oscillation persists for longer than one cycle. Simi- 
lary to Arp 102B, the estimated mass of the central black 
hole also implies a disk precession timescale of at least 
a few centuries, which is not compatible with the peri- 
ods of the _F(R)/_F(B) transient oscillations. Gilbert et 
al. (1999) found the best-fit pattern period for a spiral 
arm to produce the oscillation in i^(R)/i^(B) between 
1974 and 1998 to be greater than 25 yr. Although the 
spiral-arm model reproduces the sparsc^lv sampled trend 
of i^(R)/i^(B) to oscillate from 0.8 in 1976 to 1.5 in 1980 
and then back to ~ 0.7 after 1990, it does not match the 
profile shapes in detail. A spiral arm precessing with this 
period would predict a subsequent rise to a maximum in 
F(R)/F(B) of 1.5 in 2005, which is not consistent with 
the relatively constant F{R)/F{B) ratio of ~ 0.6 after 

2000. Thus, any regular, periodic phenomena such as 
precessing spiral arms, eccentric disks, and warps are not 
able to reproduce the irregular changes in F(R)/F(B), 
as well as the extreme changes in the shape of the wings 
of the profile and the central red component. For any of 
the above scenarios to be viable, the perturbation must 
have a life time comparable to the pattern precession. 

4.2.4. Binary Companion with Shared Accretion Disk? 

The smooth changes in the red and blue wings and 
the central component have a large amplitude, and oc- 
cur on a timescale of 5 yr. The normalized fiux in the 
blue wing increases monotonically by 50% from 1995 to 

2001, while the normalized flux in the red wing decreases 
by 65% from 1998 to 2002, and the flux in the normal- 
ized central red component decreases by 25% from 1993 
to 1998. Since a complete cycle of variation is not ob- 
served in the normalized flux of any of these regions, it is 
not yet possible to determine if these changes are in fact 
periodic or even if they repeat. If they do repeat, the re- 
currence time would be > 10 yr. A two-armed spiral does 
produce a peak in the central part of the proflle, which 
could then drift toward the blue peak (Chakrabarti & 
Wiita 1993, 1994). However, such a spiral pattern would 



12 



Gezari et al. 



also produce dramatic changes in the velocities and in- 
tensities of the red and blue peaks as the central peak 
drifts, which are not observed. The physical mechanism 
producing the observed dramatic variability must change 
the shape of the high-velocity wings of the profile, and 
the central red component, without similarly affecting 
the lower velocity gas that produces the peaks of the 
profile. One possible scenario in which changes would 
occur only in the high- velocity wings of the profile, is an 
orbiting binary companion with a common circumbinary 
accretion disk with the central black hole. The orbital 
motion of the companion perturbs the disk, and creates 
a non-axisymmetric instability that causes enhanced ac- 
cretion in the inner regions of the disk close to the com- 
panion. In addition, the binary will excite spiral waves 
in the inner parts of the circumbinary disk (Artymowicz 
& Lubow 1996; Armitage & Natarajan 2002). Both phe- 
nomena will affect the wings of the profile. For a binary 
separation less than the inner radius of the line-emitting 
portion of the disk, D < 380 r^, and a total binary mass 
of 3 X 10^ Mq, this would imply a binary separation of < 
0.005 pc, and an orbital period of < 2.2 yr. The orbital 
period of the secondary is too short to produce the > 
10 yr characteristic periods of the variability observed in 
the wings of the profile. Regardless of this, the stability 
of such a binary configuration is unlikely, since at such 
a small separation, the expected time for coalescence of 
the binary due to gravitational radiation is short, < 5 x 
10" yr (Lightman et al. 1979). 

4.2.5. Orbiting Features in the Disk? 

The lifetimes of the relatively stationary narrow fea- 
tures in the blue peak constrain whether they can be 
orbiting within the disk. The sharp spike that appears 
in 1989.50, and persists for > 2 yr at the same velocity, 
has a lifetime of at least 15% of an orbital period in the 
disk. If in orbital motion in the disk, the velocity of the 
feature would have been expected to drift by at least 2000 
km s~^ during its lifetime. The feature that appears in 
1998.27 drifts by 400 km s~^ over 1.5 months and has a 
lifetime that is short enough to be only 1% of an orbital 
period in the disk, and thus its drift of 400 km s~^ could 
be consistent with orbital motion. Another possibility 
is that the sharp features in the line profiles are shocks 
that propagate through the disk slowly (at the speed of 
sound) . Thus they can cool and dissipate quickly before 
their apparent radial velocity changes significantly. 

4.3. 3C 332 

4.3.1. Changes in the Emissivity of a Circular Disk? 

The decrease in FWHM and peak velocity separation 
in the profiles of 3C 332 during the increase in broad-line 
flux is easily accounted for by the disk model. When 
a disk is illuminated by a brighter central continuum 
source, regions farther out in the disk are photoionized, 
increasing the outer radius of the line-emitting portion 
of the disk. The increase in continuum flux can also fully 
ionize the inner regions of the disk, effectively causing an 
increase in the inner radius of the disk. The increase in 
inner and outer radius both cause lower velocity material 
to contribute more to the emission line, thus weakening 
the wings, and shifting the velocities of the peaks of the 
profile to lower velocities. A disk interpretation was also 



offered for the increase in peak velocity separation with 
decreasing broad-line flux in the double-peaked Ha pro- 
file of the LINER galaxy NGC 1097 (Storchi-Bergmann 
et al. 2003). By decreasing the inner radius of their disk 
model from 1300 to 450 r^, they could reproduce the drift 
in peaks of the Ha profile with the dimming broad-line 
flux. The decrease in inner radius was interpreted as a 
change in emissivity of the line-emitting gas, such that 
when the central continuum flux decreased, the region of 
maximum line emission shifted towards the inner regions 
of the disk. 

Figure[29lshows examples of the broad Ha profiles with 
the model disk fit for a low-state {^i ~ 190, — 540) and 
high-state (^^ = 260, = 620). Both disks have i = 36°, 
a = 1600 km s~^ , and q = 3. The spectra in the Ap- 
pendix show the least-squares scaled circular disk models 
in comparison with the profiles. Spectra before 1998 were 
fitted with the low-state model disk, and after 1998 with 
the high-state model disk. The disk profiles match well 
the decrease in FWHM and peak velocity separation, but 
the high-state profile appears to underestimate the flux 
in the very far blue wing and overestimate the flux in 
the red wing from 1998 to 2001. Beginning in 2004 there 
also appears to be an excess in the blue side of the pro- 
file, near —4000 km s~^ , that is not accounted for by 
the disk model profile. 

4.3.2. Precessmg Spiral Arm? 

The rise and faU of i^(R)/i^(B) from 1990 to 1998 
was modeled by Gilbert et al. (1999) as a one-armed 
spiral emissivity enhancement with a pattern period of 
14 yr. Although this model reproduces the behavior of 
i^(R)/i^(B) from 1990 to 1998, it does not match the 
historical i^(R)/i^(B) measurements in 1976 and 1978, 
and predicts a minimum in i^(R)/i^(B) in 2002, which 
is clearly not consistent with our continued monitoring 
of F(R)/F(B), which shows it to be nearly constant 
[F(R)/F(B) ~ 0.9] for over 8 yr after 1997 (Figure [TT|l . 
Thus the oscillation in i^(R)/i^(B) does not repeat, nor 
does it appear to be sinusoidal in shape. 

The maximum in F(R) /F{Q) in 1995 can be attributed 
to a transient excess in the red wing at -1-12,250 km 
s~^ with a FWHM ~ 3500 km s""'^ , and a maximum 
strength of 40% of the red peak. The excess in the red 
wing appears to drift to -1-10,800 km s^^ in 20 months. 
This excess may be related to the excess in the blue wing 
in 1999.46 with a width of 5000 km s~^ and a peak ve- 
locity of -11,900 km s~^ ; this could indicate a feature 
rotating with the disk that only completes one orbit be- 
fore fading away in 2002. The maximum red velocity 
of the excess of -1-12,250 km at an assumed phase 
angle of 90° and an inclination of the disk of 36°, corre- 
sponds to a bright spot at a radius of 280 . The excess 
in the blue wing appears to persist for longer than ex- 
pected for an orbiting feature, and the excess near -4000 
km s~^ does not occur when the bright spot would be 
expected to be traversing the central part of the profile. 

4.4. PKS 0235+023 
4.4.1. Circular Disk? 

The broad Ha profile of PKS 0235+023 in 1991 is 
well fitted with a circular accretion disk model with 
i = 70°, = 1020,^0 = 7000, cr = 800 km s'^ , and 
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q = 1.7 (Eracleous & Halpern 1994). The scaled disk 
model is shown in comparison with the profiles in the 
spectra in the Appendix. For all the spectra after 1991, 
the blue peak is shifted to a higher velocity, and the 
strength of the peak has increased in comparison to the 
disk model. In contrast, the red peak appears to have 
decreased in strength (and in some cases completely dis- 
appeared) , and the red wing does not extend to very high 
velocities. Even if the centroid of the entire disk profile 
is shifted blueward, the circular disk model is not a good 
fit to the profiles after 1991. 

4.5. Mkn 668 

4.5.1. Binary Broad-Line Region? 

The systematic drift of the central red peak of Mkn 
668 could be interpreted as the result of the orbital mo- 
tion of binary supermassive black holes with separate 
broad-line regions, in which the central peak is emission 
from the broad-line region of one of the orbiting black 
holes. Although the central red peak repeatedly drifts by 
-900 km s-i from 1982 to 1985, and then by -h900 km 
s^i from 1985 to 1989, and then again by > -1000 km 
s^^ from 1989 to 1995, the velocity of the peak remains 

< -1-1700 km s-^ from 1995 to 2004. The velocity curve 
does not have a regular periodic shape, and the peak ve- 
locity never passes through zero velocity, and thus in the 
case of a binary, the velocity of its center of mass would 
not be coincident with the velocity of the host galaxy. 
Thus the shift of the red peak is not consistent with or- 
bital motion of a binary broad-line region. 

4.5.2. Radiatively Accelerated Outflow? 

Marziani et al. (1993) reported a positive correlation 
between the peak velocity and the total broad-line flux 
for their spectra from 1985 to 1991. They then modeled 
the profiles as emission from a radiatively accelerated 
biconical outflow of clouds, in which the positive corre- 
lation between the peak velocity and the continuum flux 
is a consequence of radiative acceleration, and the asym- 
metry of the profile is attributed to self-absorption in the 
Balmer line or dust absorption from the backsides of the 
clouds in the outflow. Although from 1989 to 1996 the 
peak velocity and total broad-line flux decrease mono- 
tonically, we find that after 2000 this correlation breaks 
down. During a high state of broad-line flux in 2001- 
2002, during which the flux increases by a factor of 6, 
the central peak remains at a minimum peak velocity of 

< -1-1700 km s^^ . Thus the large-amplitude changes of 
the velocity of the peak are irregular, and are unrelated 
to the changes in the broad-line flux. The radiatively ac- 
celerated outflow model is thus no longer observationally 
motivated. 

4.5.3. Elliptical Accretion Disk with an Added Perturbation? 

Figure [30] shows that the shift in the red central peak 
can be modeled as a Gaussian excess at a peak veloc- 
ity of -H3300 km s^^ and a FWHM ~ 2000 km s'^ in 
1990-1995 that drifts to a peak velocity of -1800 km 
s~^ in 1996-2000, and then back to a peak velocity of 
-1-4400 km s~^ in 2003. We were able to fit the non- 
varying portion of the profile with an eccentric disk, with 
i = 35°, e< = 200, Co = 3000, g = 1.7, ti = 1500 km s"! , 



e = 0.6, and = 60° (shown in blue in Figure [30)1 . The 
maximum radial velocity of the excess and the inclina- 
tion of the disk fit can be used to determine the radius 
of the excess in gravitational radii. The maximum red- 
shifted radial velocity for a disk with an inclination of 
35° and an eccentricity of 0.6 occurs at a phase angle 
of the excess of = 25° (determined from the Doppler 
factor of an eccentric disk from equation 15 in Eracleous 
et al. 1995). Thus, for a maximum red radial velocity of 
+3300 km s"-'^ , this corresponds to a radius of 2125 r^. 
The maximum blueshifted radial velocity for an excess 
at this radius would then occur at = 255° and is -4300 
km s~^ . The observed time for the peak to shift from 
the maximum red velocity to a maximum blue velocity 
corresponds to an orbital period k, 12 yr. This model 
yields a central black hole mass from equation [4] of 1.2 x 
10* Mq. The velocity curve of a feature orbiting in the 
eccentric disk at a radius of 2125 rg is shown in compar- 
ison with the peak velocities of the Gaussian fits to the 
excess above the model disk profile in Figure [3T] The 
period of the velocity curve has been fixed to 12 yr. It is 
clear that the observed blue velocity of the feature does 
not match the amplitude of the model velocity curve, 
and when the excess returns to the red side of the disk, 
the radial velocity is too large. In general, the velocity of 
the excesses appears step-like, and not smooth, as would 
be expected for projected orbital motion. The histori- 
cal velocity curve suggests that there was an excess in 
the red side of the central peak from 1974 to 1982. This 
would be consistent with an excess that oscillates, non- 
sinusoidally, from one side of the peak to another with a 
period of ~ 12-16 yr. 

4.6. 3C 227 
4.6.1. Orbiting Bright Spot? 

The complex profile of 3G 227 cannot be modeled with 
a single disk component (circular or elliptical), since the 
profile has a low-velocity blue peak as well as a high- 
velocity blue shelf. The changes in the shape of the pro- 
file from 1990-1996 to 1999-2004 can be modeled as a 
sum of two Gaussians, one with a peak velocity of —2400 
km s^^ and a FWHM ~ 1500 km s^^ plus a Gaussian 
at a velocity of -5000 km s^^ and a FWHM ~ 2500 
km s~^ , The two Gaussians drift to the red side of the 
profile to peak velocities of -f-1900 and -f4800 km s~^ , 
respectively, from 1990 to 2000 (Figure [32| . Since the 
feature is present in the blue peak starting in 1974, this 
would imply a lifetime of the feature in the blue peak of 
> 14 yr. If the feature were orbiting from the blue side 
of the profile to the red side, then the time it takes for it 
to complete half an orbit, ss 5 yr, is inconsistent with its 
persistence in the blue peak from 1974 to 1990. Thus, 
although the feature appears to drift across the profile, 
this drift is not consistent with rotation. 

4.7. 3C 382: Elliptical Accretion Disk? 

The broad Ha profile of 3C 382 was successfully mod- 
eled with an elliptical disk by Eracleous et al. (1995), 
with i = 45°, C» = 470, Co = 3000, g = 1.7, cr = 1200 
km s^^ , e = 0.33, and (po = 110°. This disk profile 
scaled to the total broad line flux is shown in compari- 
son with the observed profiles in the spectra in the Ap- 
pendix. Although the disk model fits very well for most 
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of the observations, the drop in the red peak of the profile 
in 1994-1995, with no corresponding change in the blue 
peak, is not behavior easily produced with a disk. An or- 
biting emissivity enhancement in the disk that increases 
the strength of the red peak would be expected to tra- 
verse the profile and subsequently increase the strength 
of the blue peak. It is possible that the perturbation 
moves slowly in the azimuthal direction, and dissipates 
before it can affect the other peak. This may also be the 
case for the transient bumps that only appear in the far 
blue wing. However, the large amplitude increase in the 
red wing that subsequently occurs in the blue wing may 
be the consequence of an orbiting feature in the inner 
region of the disk, such as a bright spot or spiral arm. 
However, the excesses in the wings do not appear to tra- 
verse the lower velocity portions of the profile, which is 
not consistent with rotation. The relative constancy of 
the central part of the profile well fitted by the ellipti- 
cal disk profile can be used to place a lower limit to the 
black hole mass, since the lack of precession of the ellip- 
tical disk profile over 15 yr of the monitoring program 
implies Mbh> 1 x 10^ Mq. 

5. SUMMARY AND DISCUSSION 

The F{R)/F{B) oscillation in Arp 102B that occurred 
from 1991 to 1995, recurs in 1999 to 2004 with the same 
period. The period of the oscillations is on the dynam- 
ical timcscale, and is much shorter than the precession 
period of a spiral arm, eccentric disk, or warp. Model- 
ing the recurring oscillation as a single transient bright 
spot yields discrepant central black hole mass estimates. 
Thus, the excess in the profile that drifts from one side to 
the other may be the result of a more complex, and non- 
stable configuration of emissivity enhancements such as 
transient shocks induced by tidal perturbations. We do 
not find strong evidence for phase effects in the residuals 
of the line profiles during the oscillation of the emissivity 
enhancements from one side of the profile to the other, 
and therefore we disfavor the configuration of a collec- 
tion of clouds orbiting in a disk proposed by SOO. The 
relatively slow changes in the wings of the profile can be 
reproduced by changing the inner and outer radii of the 
line-emitting portion of the disk. In the self-consistent 
photoionized accretion disk model of Chen & Halpern 
(1989), the inner radius of the line-emitting portion of 
the accretion disk is determined by the boundary be- 
tween the thin disk and the elevated structure illuminat- 
ing the outer portions of the disk in the form of an ion 
torus or a radiatively inefficient accretion flow (RIAF; see 
Narayan et al. 1998). The outer radius is determined by 
the radius at which the disk is no longer stable to self- 
gravity, and the disk fragments. It is interesting that 
both of these radii are sensitive to the accretion rate. 
The timescale for changes in the accretion rate is on the 
viscous timescale. This timescale is orders of magnitude 
longer than the timescale of the observed variations in 
the radii of the disk. Thus changes in the accretion rate 
cannot be physically driving the changes in radii. 

The successful fit of an eccentric disk to the; non- 
variable portion of the profile of Mkn 668 makes it tempt- 
ing to attribute its oscillation of F{R) / F{B) to an orbit- 
ing emissivity enhancement. The velocity curve of the 
residuals from the disk profile model are not consistent 
with the velocity curve predicted for a feature orbiting in 



the outskirts of the disk. The historical velocity curve in- 
dicates that the oscillation of F(R)/F(B) and the central 
peak velocity is erratic, and not a regular phenomenon. 
We find that the correlation between total broad-line fiux 
and peak velocity reported by Marziani et al. (1993) is 
not compatible with the high state profiles observed in 
later spectra, during which the peak velocity remained 
at a minimum. Thus, the radiatively accelerated outflow 
model for the source of the double-peaked emission line 
is no longer observationally motivated. An orbiting bi- 
nary is also ruled out because of the irregularity of the 
peak velocity curve, and the fact that the radial velocity 
of the central peak never crosses zero velocity, requiring 
the velocity of the center of mass of the binary to not be 
coincident with that of the host nucleus. The sharp de- 
crease of the velocity of the blue peak of 3C 227 appears 
to be caused by an excess feature in the blue peak that 
traverses from one side of the profile to the other. The 
inferred orbital period of the feature of « 10 yr, however, 
is inconsistent with the persistence of the feature in the 
blue peak for over 15 yr before drifting to the red side of 
the profile. 

The decrease in FWHM and peak velocity separation 
of the profile shape of 3C 332 during its high-state of 
broad-line flux is well accounted for by the circular disk 
model, and is similar to the behavior reported in NGC 
1097, for which the peak velocity separation increased 
during a low-state in broad-line fiux. Arp 102B, on 
the other hand, demonstrates a peculiar, asymmetric re- 
sponse to flares in the broad-line flux, that are not com- 
patible with a simple response of the disk to a changing 
illumination. The repeated tendency of the blue Ha peak 
to be stronger and narrower during high-states of fiux 
also makes a random asymmetric gravitational lensing 
event an unlikely explanation. 

Sharp features that appear in the blue peak of the 
broad Ha profiles of Arp 102B and 3C 390.3 are narrow 
(< 1000 km s~^ ), and persist for months to years. The 
relatively stationary velocities of the features are incon- 
sistent with orbital motion within the disk, and require 
a well-organized velocity field. Stellar collisions with the 
accretion disk are expected to be frequent (several per 
day) and should produce a shock in the disk that can 
emit X-ray emission (Nayakshin et al. 2004), as well as 
remove a narrow plume of material. The expelled gas 
could produce the transient narrow features observed in 
the profiles, although the rate of such events is not con- 
sistent with the persistence of features over months and 
years. The appearance of the features preferentially in 
the blue peak of the profiles suggest that their blueshift 
may be a signature of an outflow. This would be expected 
from stellar collisions, since we would preferentially see 
the plumes raised over the surface of the disk. 

3C 390.3 displays a large amplitude of profile variabil- 
ity, with five coherent regions of the profile that appear 
to vary on the timescale of decades, and with no clear 
periodic correlation between each other. Even more dif- 
ficult to model is the sharpening of the blue peak, and 
the dramatic fiattening of the red peak and decrease in 
the red wing to produce a sharp shelf after 1995. Al- 
though a circular disk model provides a reasonable fit 
to the line profile in 1988, the subsequent changes in 
the profile shape make an acceptable fit to a circular 
disk model rather difficult. The excellent circular disk 
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fit to PKS 0235+023 in 1991 also breaks down for subse- 
quent spectra, in which the shape of the profile changes 
drastically. Similarly, although an elliptical disk model 
provides an excellent fit for many of the profiles of 3C 
382 during the monitoring program, the variability of 
the profile reveals inconsistencies with this simple disk 
model. The extreme variability seen in the red and blue 
wings is never detected in the lower velocity peaks of the 
profile, and although the red peak does show variability, 
the relative flux of the blue peak remains unchanged. 

6. CONCLUSIONS 

Of all the models discussed in the literature, the pho- 
toionizcd accretion disk is the most robust model for the 
origin of the double-peaked Balmer emission lines. The 
disk model provides an excellent fit to the line profiles 
with relatively few free parameters, and successfully pre- 
dicts the observed properties of the UV emission lines. 
The accretion disk model is also compatible with the spa- 
tial scales of the broad line region measured from rever- 
beration mapping experiments. 

Variability adds an extra dimension to discerning the 
structure and kinematics of the broad-line gas, and 
probes the structure of the accretion disk, the most likely 
source of the line emission. Several models for physical 
processes in the accretion disk have been proposed to ex- 
plain the ubiquitous profile variability of double-peaked 
emitters, including precessing eccentric disks and spiral 
arms, and orbiting bright spots. Individual spectra of 
the Ha profiles of six of the seven objects in this study 
are successfully fitted by line emission from a circular or 
elliptical accretion disk, however the variability of their 
line profiles reveals behavior that is not always consistent 
with the model predictions. There is no persistent, regu- 
lar variability observed in the objects that would require 
a precessing elliptical disk, one-armed spiral, or warp. 
Although the timescale of the oscillation in the relative 
flux of the peaks of the profiles in many objects appears 
to occur on the orbital timescales estimated from the 
disk model fits to the profiles, the velocity curves of the 
residuals from the model profiles are often not consis- 
tent with the projected orbital motion of a bright spot. 
In many of the objects, longer-timescale systematic vari- 
ability is seen in the wings of the profile that appears to 
be uncorrelated with the variability of the lower velocity 
regions of the profile may reflect changes in the emissiv- 
ity of the inner accretion flow. The negative correlation 
between peak velocity separation and total broad-line 
flux observed in 3C 332 is well accounted for by the disk 
model, although the asymmetric response of the profile 
during high-states of flux in Arp 102B, and possibly 3C 
390.3, is problematic. We flnd that, although the accre- 
tion disk remains the most favorable model for the origin 
of the double-peaked emission lines, most of the models 
for physical processes in the accretion disk, previously 
used in the literature to explain the variability of the 
profiles, break down when the baseline of spectroscopic 
monitoring is extended. 

Our results suggest that future efforts to model the 
variability of double-peaked emission lines should explore 
more sophisticated models. More specifically, models in 
which the perturbations have short life times, as well 



as finite onset and decay times (on the order of a few 

dynamical times, Tdyn = 6 Mg ^3^^ months) are highly 
desirable. Another family of models that warrant careful 
consideration are "stochastic models." In such models 
the observed variability is due to the collective effect of 
an ensemble of perturbations, such as shocks in the disk 
or spiral arms that are not continuous but rather are 
comprised of discrete segments. 

A promising method to probe the structure and kine- 
matics of the broad-line region is two-dimensional rever- 
beration mapping. By measuring the time delay of the 
response of the broad emission line region to the ionizing 
continuum as a function of velocity, high fidelity velocity 
delay maps can be constructed to discriminate between 
different geometries and velocity fields for the broad- 
line gas (Home et al. 2004). This method requires high 
temporal resolution, high throughput, and long and ho- 
mogenous monitoring, which would be possible with a fu- 
ture dedicated satellite mission with a UV/optical spec- 
trophotometer. In the meantime, a less demanding but 
important experiment would be to monitor a few double- 
peaked emitters on timescales shorter than the dynam- 
ical time (e.g. sampling every 1-2 weeks). Such obser- 
vations could be carried out with 2-4m ground-based 
telescopes would probe phenomena such as the devel- 
opment and propagation of shocks in the disk. We have, 
in fact, seen evidence of such phenomena in our observa- 
tions (see, e.g., ij3.2l and Figure [5]) but we were not able 
to follow them closely and study them in detail. Con- 
tinued spectroscopic monitoring of double-peaked emit- 
ters will continue to provide challenges to our current 
understanding of the nature of the broad-line region in 
double-peaked emitters and AGNs in general. 
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TABLE 1 

BLRGs IN OUR Monitoring Program 



Object 


mv 


z 


E(B-V)'' 


Arp 102B 


14.8 


0.024 


0.024 


3C 390.3 


15.4 


0.056 


0.071 


3C 332 


17.3 


0.151 


0.024 


PKS 0235+023 


17.1 


0.209 


0.033 


Mkn 668 


15.4 


0.077 


0.018 


3C 227 


16.3 


0.086 


0.026 


3C 382 


15.4 


0.059 


0.070 



^ Galactic reddening in magnitudes from Schlegel et al. (1998). 
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Fig. 1. — Temporal coverage of the spectroscopic monitoring program for the seven double-peaJjed BLRGs in our program. Each tick 
mark represents a separate observation. 
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TABLE 2 






Arp 102B: Log of Spectroscopic Observations 




UT Date 


Exp Time 


Telescope 




(sec) 




1983 Jun 06 




Palomar 5m 


1983 Jun 21 




Palomar 5m 


1985 Jun 28 




Palomar 5m 


1986 Jul 13 


3000 


Lick 3m 


1987 May 04 


1000 


Lick 3m 


1987 Aug 08 


1300 


Lick 3m 


1988 Jun 30 




Lick 3m 


1988 Jul 17 




Lick 3m 


1988 Sep 15 




Lick 3m 


1989 Apr 27 


900 


Lick 3m 


1989 May 30 




Lick 3m 


1989 Jul 01 


3600 


KPNO 2.1m 


1989 Jul 05 


835 


KPNO 2.1m 


1989 Jul 10 


900 


Lick 3m 


1990 Feb 23 


5400 


KPNO 2.1m 


1990 May 01 




Lick 3m 


1990 Jun 14 




Lick 3m 


1990 Jul 17 


1800 


Lick 3m 


1990 Jul 30 




Lick 3m 


1990 Aug 30 


1800 


Lick 3m 


1990 Nov 11 


1800 


Lick 3m 


1991 Jun 17 


9790 


KPNO 2.1m 


1991 Jun 18 


1897 


KPNO 2.1m 


1991 Jun 19 


7200 


KPNO 2.1m 


1991 Jun 20 


7200 


KPNO 2.1m 


1991 Jul 03 


7200 


KPNO 2.1m 


1991 Jul 04 


4950 


KPNO 2.1m 


1991 Aug 05 


1500 


Lick 3m 


1991 Oct 31 


1800 


Lick 3m 


1992 May 06 


7200 


KPNO 2.1m 


1992 Oct 03 


1800 


Lick 3m 


1992 Nov 19 


1200 


Lick 3m 


1993 Apr 14 


900 


Lick 3m 


1993 May 16 


3600 


KPNO 2.1m 


1993 Jun 28 


900 


Lick 3m 


1993 Jul 28 


2100 


Lick 3m 


1993 Sep 10 


3000 


Lick 3m 


1993 Sep 25 


1800 


Lick 3m 


1993 Oct 22 


1800 


Lick 3m 


1993 Nov 08 


3000 


Lick 3m 


1994 Apr 18 


1200 


Lick 3m 


1994 Jun 16 


2700 


Lick 3m 


1994 Jul 04 


3600 


KPNO 2.1m 


1994 Jul 05 


3600 


KPNO 2.1m 


1994 Jul 15 


3000 


Lick 3m 


1994 Aug 04 


3600 


Lick 3m 


1994 Sep 03 


3600 


Lick 3m 


1994 Oct 01 


4800 


Lick 3m 


1994 Nov 12 


3600 


Lick 3m 


1995 Jan 23 


1800 


KPNO 2.1m 


1995 Mar 25 


1500 


Lick 3m 


1996 Jun 02 


3600 


KPNO 2.1m 


1995 Jun 03 


3600 


KPNO 2.1m 


1995 Jun 04 


3600 


KPNO 2.1m 


1995 Sep 26 


1800 


Lick 3m 


1996 Feb 13 


3200 


KPNO 2.1m 


1996 Apr 22 


2400 


Lick 3m 


1996 Jun 13 


1711 


KPNO 2.1m 


1996 Aug 10 




Lick 3m 


1996 Sep 10 




Lick 3m 


1996 Oct 10 


2400 


Lick 3m 


1997 Feb 06 


2700 


KPNO 2.1m 


1997 Jun 08 


3600 


KPNO 2.1m 


1997 Aug 04 




Lick 3m 


1997 Sep 06 




Lick 3m 


1997 Sep 27 


3600 


KPNO 2.1m 


1998 Jan 28 


3000 


KPNO 2.1m 


1998 Apr 08 


2700 


MDM 2.4m 


1998 May 19 


3000 


KPNO 2.1m 


1998 Jun 27 


3000 


MDM 2.4m 


1998 Jul 23 




Lick 3m 


1998 Sep 21 




Lick 3m 


1998 Oct 15 


3000 


KPNO 2.1m 


1999 Jun 18 


2000 


MDM 2.4m 


1999 Sep 12 


2400 


MDM 2.4m 


2000 Jun 03 


3000 


KPNO 2.1m 


2000 Sep 25 


3600 


KPNO 2.1m 


2001 Jan 24 


1500 


MDM 2.4m 



9nni Till 9n ■^finn i^PMn 9 im 
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TABLE 3 

3C 390.3: Log of Spectroscopic Observations 



UT Date 


Exp Time 


Telescope 




(sec) 




IQS^^ Till Q 
± yoo o III y 


1200 


TCPNO dm 


1QSQ Till 01 

-Lf Of J Lll \J I. 


1800 


ivir iNv^ ^.1.111 


1 QQQ Till 09 


1 ^nn 

lOUU 


KPNO 9 1m 
IvlT^iNW ^.llU 




1800 


KPNO 9 1m 




2400 


KPNO 9 1m 

IVJT i> Zi. J-iil 


1QQ1 Tun IS 


3600 


TCPND 9 1m 

IVJT i> \J Zi. liil 


1QQ1 Tun 90 

Xyt/-L J Llli ZiL/ 


3600 


T<"PNO 9 1m 

IVJT i> Zi. liil 


1QQ1 Tun 91 

-iyt/-L J Llli Zi J. 


3600 


T<"PNO 9 1m 

IVJT i> Zi. liil 


1QQ1 Till 0*^ 
± yy -L J HI yjo 


3600 


TCPND 9 1m 

IVJT i> V_y Z.lill 


1QQ1 Till 04 


3600 


KPNO 2 1m 

xvi. 1 ^ V y z . ± ill 


1 QQ1 Till O'^ 
-Lyy± J Lll 


5000 


T<'PNn 9 1m 

IVJT iN Zi. ±111 


1 QQ9 Mav Ofi 


2700 


KPNO 9 1m 

XVJT 1^ Z[. 1.111 


1QQ4 Till (T^ 


2700 


KPNO 2 1m 


1QQ4 Tnl 0"^ 


3600 


KPNO 9 1m 

XVJT IN ^.1.111 


IQQ^ Mar 2^ 


3000 


KPNO 2 1m 


IQQ^ Tnn 04 
J-tyijO J nil U1 


3000 


KPNO 9 1m 

XVJT i'l V_/ Zr. 1.111 


IQQ^ Tnn 0"^ 

iyyo J nil UU 


3600 


KPNO 9 1m 

XVIT i\ \_/ Zr.Xlll 




ouuu 


T^PNTO 9 1m 


1 QQf{ Ar»r 91 


1200 


l_jii^lS. Olll 


IQQfi Tiin 14 
-lyyj J Llli i.'-t 


1831 


KPNO 9 1m 

IVIT i>l Zi.Xlll 


1 QQfi Or^^ 1 


1800 


T ,1 r\^' -iiT^ 
iJii^Jv Olll 


1 QQ7 FpK or 


1800 


KPNO 9 1m 

IVlT i> V_y Zi. liil 


1QQ7 Tiin OS 
xyy 1 i) Llli uo 


2700 


TCPNO 9 1m 
ivir i> Kj Zi. liil 


1 QQ7 ^on 9S 


2700 


TCPNO 9 1m 
ivir iN \j Zi. liil 


iyi70 iiJJi L/O 


1200 


MDM 2 4m 

-LVi.i_^ iV± Zj . rrlll 


1QQS Mav 20 


2540 


KPNO 2 1m 

XVi. i^ V_y Zi. J_iil 


1QQK Tun 97 
-lyyo J Llli Zi J 


1800 


MHA/T 9 4Tn 

iVll_^ivl .Z.1:lii 


1 QQ8 n*-t 1 4 


97nn 


TsTPNO 9 1m 

±\.± i > \J . ± 111 


1QQQ Tnn IS 
j-yyy j uii xo 


1800 


MDM 9 4Tn 

IVXXyiVX ^.filll 


1 QQQ Qpn 1 1 


1200 


MDM 9 4Tn 

iVXXyiVX ^.'llll 


1999 Dec 02 


1800 


KPNO 2.1m 


2000 Mar 16 


3000 


MDM 1.3m 


2000 Jun 03 


2400 


KPNO 2.1m 


2000 Sep 24 


3000 


KPNO 2.1m 


2001 Jul 23 


3600 


KPNO 2.1m 


2001 Oct 24 


2700 


KPNO 2.1m 


2002 Jun 12 


3000 


KPNO 2.1m 


2002 Oct 10 


3000 


KPNO 2.1m 


2003 Jun 08 


1800 


MDM 2.4m 


2004 Mar 28 


1800 


MDM 2.4m 


2004 Jun 10 


1800 


MDM 2.4m 


2004 Sep 21 


1800 


MDM 2.4m 
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TABLE 4 

3C 332: Log of Spectroscopic Observations 



UT Date 


Exp Time 


Telescope 




(sec) 




1 QQn A/f n^r 9S 

1 y y u J-Vi ciy Zi o 


97nn 


KPNO 9 1m 




2400 


KPNO 9 1m 

IVIT i>l V_/ .^.±111 


IQQO Tiin 19 


2700 


J-jii^iv Oill 


1 QQ 1 Foh OzL 

±TjiJ-l- -L L U Urt 


3000 


KPNO 9 1m 


1 QQ1 FpK n'^ 


3000 


KPNO 9 1m 


1QQ9 Till n? 


3200 


TCPND 9 1m 
ivir i> \j Zi. i-iii 


1QQ4 Till D'i 


3600 


KPNO 2 1m 

xvi i^ v_y Zj* j_iii 


IQQ'^ Tiin 04 


3600 


T<"PNO 9 1m 

Ivir i> V_y Z. l-iii 




3600 


TCPND 9 1m 

Ivir i> V_y Z.l-iil 


1 QQfi Anr 91 


1800 


T 1 -\ TYl 

l-jit^lv Oili 


1 QQfi Tun 1 
xyyu uu.li 


3600 


T<'PNn 9 1m 
ivir iN z.±iii 


1QQ7 Fph 07 


3600 


KPNO 2 Im 


1QQ7 Tun flQ 


3600 


KPNO 2 Im 




3600 


MDM 9 4m 

iVll_/iVX Z.^111 


1QQS Tim 97 
-Lyyo J nil ^ 1 


3600 


MDM 9 4.m 

iVlJ_/iVl Z.^111 


1 QQS Tim 98 
-Lyyo J nil 


3600 


MDM 9 4.m 

iVXl_/iVX Z.^111 


1 QQQ Tim 1 ^ 
xyyy juii ±o 


6000 


KPNO 9 1m 

XVJT i\ \_/ Z.Xlll 


9000 Tnn Ot^ 


ouuu 


T^PNTO 9 1m 

iVi i^ \_J Z. ±iil 


2000 Sep 26 


3600 


KPNO 2.1m 


2001 Jul 21 


3600 


KPNO 2.1m 


2001 Jul 22 


1008 


KPNO 2.1m 


2002 Apr 12 


3600 


MDM 2.4m 


2002 Jun 14 


3600 


KPNO 2.1m 


2002 Oct 11 


3600 


KPNO 2.1m 


2003 Jun 07 


4800 


MDM 2.4m 


2003 Jul 01 


600 


HET 9.2m 


2004 Mar 28 


6423 


MDM 2.4m 


2004 Jun 09 


7200 


MDM 2.4m 


2004 Jun 22 


600 


HET 9.2m 


2004 Aug 04 


600 


HET 9.2m 



TABLE 5 

PKS 0235+023: Log of Spectroscopic Observations 



UT Date 


Exp Time 


Telescope 




(sec) 




1991 Feb 05 


1800 


KPNO 2.1m 


1996 Oct 11 


2400 


Lick 3m 


1997 Sep 29 


3600 


KPNO 2.1m 


1998 Oct 13 


3600 


KPNO 2.1m 


1999 Sep 13 


3600 


MDM 2.4m 


1999 Dec 04 


3600 


KPNO 2.1m 


2001 Jan 24 


6000 


MDM 2.4m 


2001 Oct 25 


7200 


KPNO 2.1m 


2002 Jan 08 


3000 


MDM 2.4m 


2002 Jan 09 


3600 


MDM 2.4m 


2002 Oct 11 


3600 


KPNO 2.1m 


2002 Dec 28 


900 


HET 9.2m 


2003 Oct 23 


1800 


HET 9.2m 


2003 Dec 20 


900 


HET 9.2m 


2004 Sep 09 


1800 


HET 9.2m 


2004 Dec 13 


900 


HET 9.2m 
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TABLE 6 

Mkn 668: Log of Spectroscopic Observations 



UT Date 


Exp Time 


Telescope 




(sec) 




1QSQ Till 0*^ 
± yoy J Lii uo 


1800 


TCPNO 9 1m 


1 QQO Vnh 99 


1800 


KPND 9 1m 


1 y y u ivi dy o u 




KPNO 9 1m 


IQQ^i Tan 91 


1691 


KPNO 9 1m 

-LvJT IN V_y Zi. -LILL 


1 QQf^ ATnr 9^^ 


3600 


KPNO 9 1m 

-LVJT i> Zi. -LILL 


^QQ^ Tun OzL 


3600 


TCPND 9 1m 

-LVLT i> Zi. -LILL 


IQQfi FoV) 14 


3200 


KPNO 2 1m 


1 QQ7 PoK 07 
lyy 1 p c u u / 


9700 


T^PNO 9 1m 
1 vi IN W 1111 


1QQ7 Tun OQ 
-Lyy 1 u uLi V-jy 


2700 


KPNO 9 1m 

-LVLT i> \J Zi.±lLL 




3600 


KPNO 2 1m 

XVX 1 ^ V_7 Zj • -L ILl 


1998 Apr 08 


1800 


MDM 2.4m 


1998 May 19 


2700 


KPNO 2.1m 


1998 Jun 27 


2000 


MDM 2.4m 


1999 Jun 15 


2700 


KPNO 2.1m 


2000 Jun 03 


2000 


KPNO 2.1m 


2001 Jan 24 


1800 


MDM 2.4m 


2001 Jul 21 


3000 


KPNO 2.1m 


2002 Jan 08 


1800 


MDM 2.4m 


2003 Jan 26 


4400 


MDM 1.3m 


2004 Mar 29 


1800 


MDM 2.4m 


2004 Jun 09 


1800 


MDM 2.4m 





TABLE 7 






3C 227: Log of Spectroscopic Observations 




UT Dato 


Ex]) Time 


Telescope 




(«■'■) 




1990 Feb 22 


7200 


KPNO 2.1m 


1993 Dec 13 


1500 


KPNO 2.1m 


1995 Jan 23 


3600 


KPNO 2.1m 


1995 Mar 24 


3600 


KPNO 2.1m 


1996 Feb 15 


3200 


KPNO 2.1m 


1997 Feb 06 


3600 


KPNO 2.1m 


1998 Jan 29 


3600 


KPNO 2.1m 


1998 May 20 


2700 


KPNO 2.1m 


1998 Dec 20 


2400 


MDM 2.4m 


1999 Dec 04 


3600 


KPNO 2.1m 


2001 Jan 24 


4800 


MDM 2.4m 


2002 Jan 08 


2700 


MDM 2.4m 


2002 Jan 11 


2700 


MDM 2.4m 


2002 Apr 12 


2400 


MDM 2.4m 


2002 Dec 27 


300 


HET 9.2m 


2003 Jan 26 


4800 


MDM 1.3m 


2003 Apr 05 


300 


HET 9.2m 


2003 Oct 31 


300 


HET 9.2m 


2003 Dec 22 


300 


HET 9.2m 


2004 Mar 28 


2700 


MDM 2.4m 



22 



Gezari et al. 



TABLE 8 

3C 382: Log of Spectroscopic Observations 



UT Date 


Exp Time 


Telescope 




(sec) 




1 QQQ Tnl 09 


1 ^nn 


T*rPNO 9 1m 


1QSQ Till OzL 


1800 


KVNCl 9 1m 


1 QQQ Till r\r^ 


1 9nn 


T^TPNO 9 1m 


1 QQO Mp\\r '\f\ 

-LyfJ ivicty rj\j 


2400 


KPNO 9 1m 

IVJT IN V_y Zi. llll 


1QQ1 Tun 18 

o nil _10 


3600 


KPNO 2 1m 

xvi IN \_y Zj* j_iii 


Xi/iJ-L J nil ZiL/ 


3600 


TCPND 9 1m 

IVJT 1> Zi. llll 


1QQ1 Tun 91 

Xi/t/J- J Llll Zi J. 


3600 


T<"PNO 9 1m 

IVJT 1> v_y Zi. llll 


1 QQ/I Till 

iyy4: J Lll uo 


91 on 


T^PNO 9 1m 
1 vi IN W llll 


^QQ^ Tun ^'^ 
lyyu J nil uo 


1800 


T^PNO 9 1m 

IVJT 1> v_y Zi. llll 


IQOfi Anr 21 

± T / ( ./l. L/l ^ ±. 


1200 


T .If li" "-im 

iJl^'XV ijlll 


1 QQf^ Anr 91 


1200 


T .1 r'hr -im 
J_jHjJv Olll 


IQQfi TiiTi 14 


442 


KPNO 2 1m 

x\.± 1 ^ Z/ ■ xiii 


IQQfi Ort in 


900 


T .1 r*!;" "^m 

UlL'lV <_llll 


1 QQ7 Mar 91 


600 


MDM 9 4Tn 


1 QQ7 Tnn flR 

-Li7£7 1 J nil UO 


1200 


KPNO 9 1m 

XVJT Zi.Xlll 


1 QQ7 CJpT-i 9Q 


1800 


KPNO 9 1m 

XVJT i'l Zr.Xlll 


1 QQa A r>r OS 


900 


1VXX_/1VX ^.^111 


1 QQQ Tnn 97 

lyyo J Llll ^1 


1 9nn 


iVl LJ IVl Z . Ill 1 1 


1 QQ^^ Ciri 1 
lyyo v^L-L ±u 


1800 


KPNO 9 1m 

XVX In Zr.Xlll 


1 QQQ Tun 1 


1200 


KPNO 2 1m 

XVX IN V_y Zf . X 111 


1 QQQ ^on 1 1 

lyfy o*--]-j IX 


900 


IVXXyiVX Zi.^lll 


1 QQQ T)or 0*^ 
lyyy i-/i^t, uo 


1800 


KPNO 9 1m 

XVX IN KJ Zi. Xlil 


ofinn Tiin nzi 

ZUUU J Llll U^ 


1 ^nn 


T^PNTO 9 1m 

XVX IN kJ Z. llll 


9000 ^on 9^^ 

^UUU 0(-]J 


1800 


TCPNO 9 1m 
XvX IN V_y Zi. Xlil 


2001 Jul 21 


1800 


KPNO 2.1m 


2001 Oct 23 


1200 


KPNO 2.1m 


2002 Apr 12 


1200 


MDM 2.4m 


2002 Jun 03 


300 


HEX 9.2m 


2002 Jun 13 


2100 


KPNO 2.1m 


2002 Oct 10 


1800 


KPNO 2.1m 


2003 Apr 01 


300 


HET 9.2m 


2003 Jun 07 


3600 


MDM 2.4m 


2004 Mar 29 


1200 


MDM 2.4m 


2004 Jun 09 


1200 


MDM 2.4m 


2004 Jun 25 


300 


HET 9.2m 


2004 Sep 09 


150 


HET 9.2m 
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TABLE 9 
Arp 102B: Broad Hq Properties 



UT Date 


Flux 


V(B) 


P(B) 


F(R)/F(B) 


FWQM 


QMS 


QMC 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1983 Jun 06 


21.0±0.6 


-5.0 


0.069±0.002 


0.98 


16.9 


0.9 


0.9 


1983 Jun 21 


21.0±0.6 


-5.0 


0.069±0.002 


0.98 


16.9 


0.9 


0.9 


1985 Jun 28 


16.6±0.6 


-4.8 


0.059±0.002 


0.80 


17.2 


0.7 


0.4 


1986 Jul 13 


19.1±0.8 


-5.2 


0.079±0.003 


0.72 


16.4 


0.9 


0.4 


1987 May 04 


23±1 


-4.8 


0.087±0.005 


0.76 


16.5 


0.9 


0.4 


1987 Aug 08 


20.3±0.8 


-4.9 


0.077±0.003 


0.72 


16.8 


1.0 


0.4 


1988 Jun 30 


13.8±0.7 


-5.4 


0.060±0.003 


0.70 


16.3 


0.8 


0.1 


1988 Jul 17 


13.9±0.9 


-5.2 


0.056±0.004 


0.70 


16.6 


0.7 


0.1 


1988 Sep 15 


16.0±0.7 


-5.3 


0.065±0.003 


0.71 


16.6 


0.7 


0.1 


1989 Apr 27 


17.1±0.9 


-5.0 


0.068±0.004 


0.74 


17.2 


1.2 


0.7 


1989 May 30 


14.7±0.7 


-5.3 


0.057±0.003 


0.74 


18.0 


1.2 


0.7 


1989 Jul 01 


14.4±0.7 


-5.2 


0.056±0.003 


0.74 


17.8 


1.3 


0.7 


1989 Jul 05 


19.8±0.6 


-5.2 


0.085±0.003 


0.68 


16.3 


0.8 


0.2 


1989 Jul 10 


18±2 


-5.2 


0.068±0.007 


0.73 


18.3 


1.1 


0.6 


1990 Feb 23 


14.8±0.6 


-5.3 


0.052±0.002 


0.70 


18.3 


1.0 


0.5 


1990 May 01 


13.7±0.6 


-5.0 


0.045±0.002 


0.72 


19.0 


0.8 


0.3 


1990 May 30 


15.1±0.9 


-5.0 


0.051±0.003 


0.70 


19.0 


0.8 


0.3 


1990 Jun 14 


13.8±0.8 


-5.3 


0.049±0.003 


0.64 


18.8 


0.7 


0.1 


1990 Jul 17 


17±1 


-5.1 


0.054±0.003 


0.74 


18.7 


1.2 


0.8 


1990 Jul 30 


12.4±0.7 


-5.1 


0.044±0.003 


0.65 


18.5 


1.1 


0.4 


1990 Aug 30 


14.7±0.7 


-5.0 


0.050±0.002 


0.75 


18.3 


1.0 


0.6 


1990 Nov 11 


16±1 


-4.7 


0.053±0.003 


0.84 


17.6 


1.0 


0.8 


1991 Jun 17 


12.5±0.5 


-5.4 


0.050±0.002 


0.62 


17.9 


0.7 


0.0 


1991 Jun 18 


13.6±0.5 


-5.5 


0.053±0.002 


0.64 


18.2 


0.9 


0.1 


1991 Jun 19 


13.8±0.5 


-5.5 


0.053±0.002 


0.64 


18.2 


0.9 


0.1 


1991 Jun 20 


14.6±0.6 


-5.5 


0.057±0.002 


0.60 


18.3 


0.7 


-0.1 


1991 Jul 03 


12.9±0.5 


-5.5 


0.051±0.002 


0.64 


18.0 


0.7 


0.0 


1991 Jul 04 


12.4±0.5 


-5.5 


0.050±0.002 


0.61 


18.2 


0.8 


0.0 


1991 Aug 05 


16.9±0.5 


-5.2 


0.062±0.002 


0.66 


18.6 


1.0 


0.3 


1991 Oct 31 


12.9±0.7 


-5.4 


0.047±0.002 


0.63 


17.7 


0.4 


-0.2 


1992 May 06 


17.0±0.8 


-5.2 


0.056±0.003 


0.85 


19.0 


1.4 


0.9 


1992 Oct 03 


13.0±0.6 


-5.2 


0.048±0.002 


0.83 


17.9 


1.0 


0.7 


1992 Nov 19 


15±1 


-5.1 


0.048±0.003 


0.84 


18.6 


1.0 


0.7 


1993 Apr 14 


14.0±0.7 


-5.1 


0.054±0.003 


0.63 


17.2 


0.9 


0.1 


1993 May 16 


13.9±0.8 


-5.0 


0.051±0.003 


0.69 


17.3 


1.1 


0.4 


1993 Jun 28 


16.4±0.6 


-5.0 


0.059±0.002 


0.74 


17.4 


0.9 


0.4 


1993 Jul 28 


14.9±0.7 


-5.1 


0.059±0.003 


0.65 


17.4 


0.9 


0.2 


1993 Sep 10 


15.7±0.6 


-5.2 


0.062±0.002 


0.64 


17.4 


0.8 


0.1 


1993 Sep 25 


18.2±0.8 


-5.1 


0.073±0.003 


0.66 


17.4 


0.8 


0.1 


1993 Oct 22 


15.1±0.6 


-5.2 


0.062±0.002 


0.67 


17.5 


0.7 


0.3 


1993 Nov 08 


15.4±0.6 


-5.3 


0.060±0.002 


0.65 


17.7 


1.0 


0.3 


1994 Apr 18 


14.1±0.6 


-5.1 


0.056±0.002 


0.76 


17.4 


1.2 


0.9 


1994 Jun 16 


17.3±0.7 


-5.3 


0.068±0.003 


0.79 


17.6 


1.3 


0.9 


1994 Jul 04 


15.1±0.6 


-5.4 


0.059±0.002 


0.80 


17.4 


1.2 


0.9 


1994 Jul 05 


15.8±0.6 


-5.4 


0.059±0.002 


0.85 


17.4 


1.1 


0.9 


1994 Jul 15 


16.2±0.8 


-5.3 


0.061±0.003 


0.83 


17.5 


1.2 


0.9 


1994 Aug 04 


18.7±0.7 


-5.3 


0.068±0.003 


0.86 


17.3 


1.0 


0.9 


1994 Sep 03 


15.9±0.8 


-5.3 


0.062±0.003 


0.78 


17.1 


1.0 


0.7 


1994 Sep 08 


15.9±0.8 


-5.3 


0.062±0.003 


0.78 


17.1 


1.0 


0.7 


1994 Oct 01 


15.3±0.7 


-5.3 


0.058±0.002 


0.82 


17.4 


1.1 


0.8 


1994 Nov 12 


14.5±0.9 


-5.3 


0.058±0.004 


0.71 


17.3 


0.9 


0.5 


1995 Jan 23 


14.8±0.6 


-5.2 


0.056±0.002 


0.72 


16.9 


0.8 


0.5 


1995 Mar 25 


15.3±0.7 


-5.2 


0.056±0.003 


0.82 


17.4 


1.2 


0.9 


1995 Jun 02 


15.2±0.6 


-5.1 


0.055±0.002 


0.71 


17.6 


1.0 


0.5 


1995 Jun 03 


16.5±0.5 


-5.1 


0.057±0.002 


0.75 


17.7 


1.1 


0.7 


1995 Jun 04 


16.5±0.5 


-5.1 


0.057±0.002 


0.75 


17.7 


1.1 


0.7 


1995 Jul 06 


13.4±0.8 


-5.1 


0.048±0.003 


0.72 


18.0 


1.1 


0.4 


1995 Sep 26 


12.7±0.7 


-5.0 


0.051±0.003 


0.64 


17.6 


1.2 


0.3 


1996 Feb 13 


16.3±0.6 


-4.8 


0.060±0.002 


0.70 


17.4 


0.7 


0.3 


1996 Apr 22 


12.8±0.3 


-4.8 


0.051±0.001 


0.64 


17.7 


1.2 


0.4 


1996 Jun 13 


14.2±0.8 


-5.1 


0.053±0.003 


0.74 


17.8 


1.3 


0.7 


1996 Aug 10 


17±1 


-5.2 


0.066±0.004 


0.69 


17.4 


1.1 


0.4 


1996 Sep 10 


13.4±0.6 


-5.3 


0.059±0.003 


0.62 


17.4 


1.0 


0.2 


1996 Oct 10 


14.3±0.7 


-5.3 


0.060±0.003 


0.64 


16.9 


0.9 


0.3 


1997 Feb 06 


12.7±0.6 


-5.3 


0.056±0.003 


0.65 


17.1 


0.5 


-0.1 


1997 Jun 08 


14.1±0.6 


-5.2 


0.055±0.002 


0.69 


17.9 


0.8 


0.3 


1997 Aug 04 


15.2±0.8 


-5.3 


0.062±0.003 


0.65 


17.2 


0.5 


0.0 


1997 Sep 06 


16±1 


-5.3 


0.066±0.005 


0.64 


16.6 


0.6 


-0.1 


1997 Sep 27 


18.5±0.7 


-5.3 


0.070±0.003 


0.78 


16.7 


0.7 


0.3 


1998 Jan 28 


21±1 


-5.0 


0.079±0.005 


0.85 


16.7 


0.9 


0.6 


1998 Apr 08 


21.0±0.9 


-5.1 


0.083±0.004 


0.73 


16.4 


0.7 


0.2 


1998 May 19 


19.8±0.6 


-5.2 


0.085±0.003 


0.68 


16.3 


0.8 


0.2 


1998 Jun 27 


31±2 


-5.2 


0.141±0.007 


0.70 


15.9 


0.7 


0.2 


1998 Jul 23 


23±2 


-4.8 


0.098dz0.006 


0.67 


16.0 


0.7 


0.1 


1998 Sep 21 


17.7±0.8 


-4.9 


0.078±0.004 


0.58 


15.8 


0.5 


-0.4 


1998 Oct 15 


23.3±0.8 


-4.9 


0.092±0.003 


0.68 


16.4 


0.7 


0.1 


1999 Jun 18 


17.0±0.7 


-5.1 


0.075±0.003 


0.67 


15.7 


0.4 


-0.1 


1999 Sep 12 


17.0±0.6 


-5.0 


0.077±0.003 


0.62 


15.9 


0.4 


-0.3 


9nnn Tun ns 


1 fi s-t-n 7 


n 


n nfi7-i-n nn-i 


n 77 


Ifi 8 


1 n 


n 
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TABLE 10 
3C 390.3: Broad Ha Properties 



UT Date 


Flux 


V(R) 


V(B) 


P(R) 


P(B) 


F(R)/F(B) 


FWQM 


QMS 


QMC 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1 QQQ Till Q 


9Q-I-'^ 
zy^o 


4 




n OQ-I-O 01 


n 1 9-1-0 01 

U. -LZ3IU.U-L 


n 77 


1 Q 


8 
U.o 


^ 
u.o 


1QSQ Till 01 

j-ijOU Lll \J 1. 


45ibl 


4.1 


—3.6 


n 1 'iR-l-O OO'i 


n 0074-0 00"^ 


0.63 


13.5 


0.4 


— 0.1 


1QSQ Till 09 


50ibl 


4.1 


—3.5 


n 1 p;o-ho OO'i 


99A-kO 00^^ 


0.63 


13.3 


0.3 


—0.1 


1 QQO T^oK 99 


oo^^ 


4 2 


4 


n 1 R-i-0 01 


n 9Q-U0 09 
u. zozcu.uz 


u.uo 


lO.O 


1 
U. J- 


9 
u. z 


1 QQO Mp\\r '^0 


45ibl 


4.1 


—3.2 


fi 1 QA-kO OO'^ 


911 -UO OO'^ 
U. Z ± 1- ^u . uuu 


0.59 


13.0 


0.2 


—0.2 


1QQ1 Tun IS 


48ib2 


4.3 


—3.1 


n 1 47-kO 007 
u. i-i 1 inu.uu / 


n 91 -UO 01 

U.Z 1. ^U . U 1. 


0.68 


13.1 


0.1 


—0.1 


1QQ1 Tun 9n 


47ibl 


4.4 


—3.1 


n 1 4*^-^0 OO'^ 
u. iioinu . uuo 


90S-U0 OO'^ 


0.68 


13.4 


0.3 


0.0 


1QQ1 Tun 21 


51ib2 


4.3 


—3.2 


1 ^9+n 007 


99-hO 01 


0.67 


13.5 


0.4 


0.0 


1QQ1 Till 0*^ 

±tyiJ 1. O Lll W O 


50ib2 


4.4 


—3.2 


n 1 c;fi-kO 007 


n 99-UO 01 

U.ZZ^U.U-L 


0.70 


13.4 


0.5 


0.1 


1QQ1 Till 04 
±yy± u U-i 


50di3 


4.3 


—3.2 


n 1 4S-I-0 OOS 


n 99-Un 01 

U.ZZ^U.U-L 


0.67 


13.3 


0.3 


0.0 


1QQ1 Till 0^ 


52di3 


4.4 


—3.2 


*J. 1 t n t 1 1 I.I n 1(1 


9^-1-0 01 


0.68 


13.5 


0.3 


0.0 


1 QQ9 Mav nfi 


58ib5 


4.2 


—3.4 


n 1 e-un 01 


n 9f;-|-n 09 


0.68 


13.0 


0.1 


—0.2 


1QQ4 Till 0*^ 


35-1-2 


4.2 


—3.3 


nil c_|_r) nny 


n 1 7-i-n ni 


0.68 


13.4 


0.5 


0.0 


1QQ4. Till Cl^ 


34ib2 


4.2 


—3.4 


nil Q-Ln 00^ 


n 1 K7-|-n 007 

U. -LU 1 ^U.UU 1 


0.69 


13.4 


0.5 


0.1 


1 QQi^; Mar 9^ 


42ib3 


4.4 


—3.3 


141 -1-0 OOQ 


n iQ-i-0 01 


0.74 


14.5 


1.1 


0.5 


IQQ'S Inn 04 

J^iytyvj tJLi.ll \Jrr 


38di2 


4.3 


—3.4 


n 1 s4-i-n nn^ 


n 1 «n-|-n 007 


0.74 


14.1 


0.9 


0.3 


IQQf^ Tun 0'^ 

j.yyo J U.11 uu 


42ib2 


4.2 


—3.3 


n iae;4-0 007 


n 1 Qfi-l-O OOQ 


0.76 


14.3 


1.0 


0.4 


1 QQfi FpVi 1 


56ib2 


4.0 


—2.7 


n 1 Q^-i-0 008 


91 S-UO OOQ 
u. z j_ozizu. uuy 


0.82 


14.2 


0.6 


0.3 


1 QQ« Anr 91 


47-H 


4 n 


9 Q 


n 1 1^1 -Ufi 004 


1 Q94-0 00^^ 
u. lyznzu.uuo 


n 78 


lO.U 




u.o 




u.u 


1 QQiS Or^^ 1 


49ib3 


3.8 


— 2.9 


n 1 j^-un 01 


n 9q-U0 01 

U . ZO^U. U 1. 


0.69 


13.6 


0.4 


0.0 


1 QQ7 Foh OR 


47ib4 


3.3 


— 3.0 


fl 1 4-1-0 01 


n 99-UO 09 

U. ZZZIZU. uz 


0.61 


13.3 


0.3 


—0.2 


1QQ7 Tiin OS 
± yy / J Llli uo 


40ib2 


3.3 


—3.2 


n 1 1 7-Un 00*^ 


n 1 Q-UO 01 
u . -Ly mu . U-L 


0.60 


14.3 


0.7 


0.0 


1 QQ7 ^on 9K 
±yy/ ocjj zo 


41ibl 


3.2 


—3.2 


191 4-0 00'^ 

VJ. 1. 1. _I_U . lJUO 


1 Q'^-UO 00"^ 
u. lyomu.uuo 


0.60 


14.4 


0.8 


0.0 


1 QQK Anr flK 

±yyo jT-jJi uo 


43dil 


3.6 


—3.2 


n 1 94-1-0 00'^ 


1 QS-UO 00"^ 
u. lyo^u.uuo 


0.64 


14.4 


0.6 


—0.1 


1 QQS A/Tav 90 
1- y y o ivi cty z u 


49i6 


3.6 


—3.1 


n 1 f;4-o 09 

U. 1- o_izu. uz 


A 99-UO 0*^ 
u.zz^u. uo 


0.66 


14.7 


0.8 


0.0 


1QQS Tun 27 

1-^I/^yO U Llli ^ 1 


58ib5 


3.8 


—3.2 


n 1 S-l-0 01 


9^-hO 09 

U. iiU_l_U . Uii 


0.71 


14.9 


0.8 


0.2 


1QQK Ori 14 
lyyo v_/(_-Lr ±'-t 


51ib4 


3.2 


—3.4 


n 1 4-1-0 01 


n 944-0 09 
u . iiTi j[_u . uz 


0.58 


14.7 


0.6 


—0.2 


-Li/C/t/ U Llli iO 


40di3 


2.8 


—3.3 


n 1 1 fi-l-0 OOQ 


n 20+0 02 


0.53 


14.9 


0.7 


—0.2 


1 QQQ QpT-i 1 1 


40zb2 


3.0 


—3.4 


Oil 0-1-0 OOR 


n iQ-|-n 01 


0.54 


14.9 


0.5 


—0 


1999 Dec 02 


39±4 


3.0 


-3.4 


0.12±0.01 


0.19±0.02 


0.59 


14.9 


0.4 


-0.3 


2000 Mar 16 


46±4 


3.7 


-3.4 


0.14±0.01 


0.22±0.02 


0.65 


15.1 


0.4 


-0.2 


2000 Jun 03 


49±4 


3.6 


-3.5 


0.14±0.01 


0.23±0.02 


0.62 


14.8 


0.3 


-0.3 


2000 Sep 24 


49±2 


3.7 


-3.6 


0.139±0.007 


0.24±0.01 


0.61 


14.0 


0.0 


-0.4 


2001 Jul 23 


42±3 


3.7 


-3.5 


0.127±0.008 


0.21±0.01 


0.64 


14.8 


0.4 


-0.1 


2001 Oct 24 


34±2 


3.0 


-3.5 


0.113±0.006 


0.19±0.01 


0.61 


13.9 


-0.3 


-0.7 


2002 Jun 12 


32±2 


3.4 


-3.1 


0.108±0.006 


0.157±0.009 


0.69 


14.3 


-0.1 


-0.3 


2002 Oct 10 


41=b2 


3.5 


-3.3 


0.142±0.008 


0.21±0.01 


0.70 


13.5 


-0.1 


-0.3 


2003 Jun 08 


44±4 


3.5 


-3.2 


0.14±0.01 


0.20±0.02 


0.66 


13.8 


-0.1 


-0.5 


2004 Mar 28 


66±8 


3.4 


-3.7 


0.19±0.02 


0.36±0.04 


0.56 


13.3 


-0.1 


-0.7 


2004 Jun 10 


57i:4 


3.7 


-3.3 


0.18±0.01 


0.30±0.02 


0.61 


13.7 


0.2 


-0.4 


2004 Sep 21 


63±5 


3.6 


-3.3 


0.19±0.01 


0.33±0.03 


0.61 


13.8 


0.1 


-0.5 



Note. — Col. (1): Date of observation. Col. (2): Flux of broad Ho relative to [O I]A6300 with la errors. Col. (3-4); Flux-weighted 
velocity centroid of top 10% of the red peak and blue peak in units of 10'^ km . Col. (5-6); Average flux in top 10% of the red peak and 
blue peak relative to [O I] A6300 with la errors. Col. (7): Ratio of average flux in fixed velocity intervals corresponding to the top 10% of the 
red peaJj and blue peak in the mean profile. Col. (8): Pull width at quarter- maximum in units of 10^ km s~^. Col. (9): Quarter-maximum 
shift measured from central velocity at quarter maximum in units of 10^ km s~^. Col. (10): Flux-weighted centroid at quarter-maximum 
in units of 10^ km 
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TABLE 11 
3C 332: Broad Ha Properties 



UT Date 


Flux 


V{R) 


V(B) 


P{R) 


P(B) 


F{R)/F(B) 


FWHM 


HMS 


KMC 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1 QQn A/f 9S 
1 y y u ivi cny a o 




o.y 


7 A 


n 090-1-0 noi 


094 +0 001 


U.U4 


21 8 


1 

1.0 


u.u 


J- y y J. VI t%y O y) 


9.7i0.6 


11.4 


—7.5 


n 097-1-0 009 


n o'^'^+O 009 


0.66 


22.4 


1.3 


0.4 


IQQO Tun 19 


12.1ib0.4 


9.4 


—8.3 


O'^'^-l-O 001 


04S+0 009 


0.68 


20.2 


9 


0.2 


1 QQ 1 PoK OA 


^ 7-1-0 


o.u 




01 9-1-0 001 


01 7+0 001 


u.ou 


9"^ 
zo.o 


1 


]^ ]^ 


1 QQ1 FpK (T^ 


5 9-j-Q 3 


9.6 


— 7.6 


01 9-i-O 001 


n 01 S+0 001 

U. U -LO_L_U. UU -L 


0.72 


22.6 


1.5 


0.8 


1QQ2 Till (17 


6.6±0.4 


10.2 


— 7.5 


n 01 4+0 Oni 


ni 7+0 001 


0.88 


23.9 


1.8 


1.6 


1QQ4 Till (T^ 


6.4it0.4 


9.9 


— 7.0 


01 S-i-0 001 


01 '^+0 001 


1.20 


23.7 


1.8 


2.6 




10.8it0.4 


11.2 


—6.0 


09S+0 oni 


097+0 001 


1.01 


24.2 


2.1 


2.5 


1 QQfi FpK 1 R 


ll.litO.4 


10.1 


—6.9 


097-1-0 001 


09S+0 001 

U. UZrO^U.UUX 


0.90 


24.0 


1.6 


1.8 




9.7it0.4 


10.2 


—6.9 


09R-I-0 001 


09^^+0 001 
u.u.^u^u.uux 


0.86 


24.0 


1.5 


1.5 


IQQfi Inn 1 'i 


10.2±0.4 


9.8 


—7.1 




094+0 001 


1.01 


24.5 


1.6 


2.0 


1 QQ7 TTfsh 07 


10.2±0.4 


10.0 


—7.4 




n 097+0 noi 


0.84 


23.3 


1.5 


1.4 


1QQ7 Tun flQ 


10.4±0.4 


9.6 


—7.1 


n 024+0 noT 


n 028+0 noi 

U . U tiO _1_ U • UU X 


0.80 


22.8 


1.2 


1.0 




17.0±0.6 


8.0 


—7.1 


n n^fi+o noi 


n o4fi+o no9 


0.81 


21.6 


0.9 


0.5 


1QQS Tnn 97 
-Lyyo J nil ^ I 


18.5it0.7 


7.9 


—7.1 


n oAl +0 009 


n OAQ+n n09 


0.83 


21.7 


0.8 


0.4 


1QQ« TiiTi 98 

-Li7t70 U Ull i;0 


16.1±0.6 


8.2 


—7.0 


n nsfi+0 noi 


041 _i_n 002 

U.U^X_l_U.UUii 


0.87 


21.9 


0.8 


0.5 


1 QQQ Tnn 1 ^ 
xyyy juii xu 


36ib2 


8.1 


—6.0 


o^^'^+o 00^ 


n n«Q+n nOR 

U . UOO^U . UUO 


0.90 


20.6 


0.5 


0.3 


9nnn Tun 


32-j-i 


7.4 


—5.5 


079+0 00'^ 


0^^0+0 OO'i 

U.UOUZCU.UUO 


0.88 


20.8 


0.1 


0.4 


2000 Sep 26 


33.4±0.6 


7.5 


-5.6 


0.074±0.001 


0.086±0.002 


0.85 


20.0 


0.6 


0.3 


2001 Jul 21 


18.3±0.6 


8.8 


-5.6 


0.041±0.001 


0.049±0.002 


0.77 


20.2 


0.4 


0.0 


2001 Jul 22 


26.4±0.8 


8.0 


-5.3 


0.061±0.002 


0.077±0.002 


0.76 


19.7 


0.6 


0.1 


2002 Apr 12 


33±3 


7.5 


-5.4 


0.076±0.006 


0.076±0.006 


0.93 


20.0 


0.6 


0.2 


2002 Jun 14 


18.9±0.6 


7.9 


-5.8 


0.041±0.001 


0.045±0.002 


0.88 


20.7 


0.6 


0.3 


2002 Oct 11 


26±1 


7.9 


-6.2 


0.059±0.002 


0.063±0.002 


0.89 


20.8 


0.7 


0.4 


2003 Jun 07 


43±4 


7.3 


-5.7 


O.lliO.Ol 


0.12±0.01 


0.81 


21.1 


0.2 


0.1 


2003 Jul 01 


49±2 


7.3 


-5.3 


0.120±0.005 


0.115±0.004 


0.95 


19.8 


0.8 


0.5 


2004 Mar 28 


35±3 


7.1 


-5.0 


0.072±0.005 


0.095±0.007 


0.79 


20.0 


1.2 


0.6 


2004 Jun 09 


59±6 


7.6 


-5.2 


0.14±0.01 


0.16±0.02 


0.81 


19.6 


1.0 


0.3 


2004 Jun 22 


46±2 


7.6 


-5.0 


0.106±0.004 


0.119±0.005 


0.91 


20.0 


1.1 


0.7 


2004 Aug 04 


37±1 


7.9 


-5.1 


0.083±0.003 


0.096±0.003 


0.87 


19.9 


1.2 


0.8 



Note. — Col. (1): Date of observation. Col. (2): Flux of broad Ha relative to narrow Hq with la errors. Col. (3-4): Flux-weighted 
veloeity centroid of top 10% of the red peak and blue peak in units of 10'^ km s^^. Col. (5-6): Average flux in top 10% of the red peak and 
blue peak relative to narrow Ha with Icr errors. Col. (7): Ratio of average flux in fixed veloeity intervals eorresponding to the top 10% of 
the red peak and blue peak in the mean profile. Col. (8): Full width at half-maximum in units of 10'^ km s^-"^. Col. (9): Half-maximum 
shift measured from central velocity at half maximum in units of 10"^ km s~^. Col. (10): Flux- weighted centroid at half-maximum in units 
of 103 km s-l 



TABLE 12 
PKS 0235-1-023: Broad Ha Properties 



UT Date 


Flux 


V(B) 


P(B) 


F(R)/F(B) 


FWQM 


QMS 


QMC 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1991 Feb 05 


9.0±0.9 


-3.6 


0.034±0.003 


1.37 


14.7 


0.3 


0.2 


1996 Oct 11 


7.9±0.8 


-4.3 


0.030±0.003 


0.70 


14.6 


-1.3 


-1.3 


1997 Sep 29 


3.8±0.4 


-4.4 


0.020±0.002 


0.52 


13.5 


-1.5 


-2.1 


1998 Oct 13 


7.1±0.7 


-4.7 


0.042±0.004 


0.52 


12.0 


-1.0 


-1.7 


1999 Sep 13 


12±1 


-4.7 


0.059±0.006 


0.56 


12.9 


-1.0 


-1.5 


1999 Dec 04 


8.4±0.8 


-4.5 


0.046±0.004 


0.67 


13.5 


0.1 


-0.9 


2001 Jan 24 


7.3±0.7 


-4.1 


0.036±0.003 


0.67 


13.4 


-0.8 


-1.4 


2001 Oct 25 


15±1 


-4.7 


0.066±0.006 


0.71 


13.8 


-0.7 


-1.2 


2002 Jan 09 


11±1 


-4.6 


0.052±0.005 


0.70 


14.5 


-0.4 


-1.1 


2002 Oct 11 


17±2 


-4.5 


0.066it0.006 


0.77 


14.7 


-0.7 


-1.0 


2002 Dec 28 


8.2±0.8 


-4.2 


0.036±0.003 


0.75 


14.4 


-1.0 


-1.4 


2003 Oct 23 


13±1 


-4.6 


0.051±0.005 


0.76 


15.7 


-1.5 


-1.5 


2003 Dec 20 


12±1 


-4.9 


0.046±0.004 


0.74 


15.2 


-1.6 


-1.5 


2004 Sep 09 


12±1 


-4.9 


0.054±0.005 


0.72 


14.3 


-1.1 


-1.4 


2004 Dec 13 


5.4±0.5 


-5.0 


0.028±0.003 


0.68 


13.8 


-1.3 


-1.7 



Note. — Col. (1): Date of observation. Col. (2): Flux of broad Ho relative to narrow Ha with Icr errors. Col. (3): Flux-weighted 
velocity centroid of top 10% of blue peak in units of 10^ km s~^. Col. (4): Average flux in top 10% of blue peak relative to narrow Ha 
with 1(T errors. Col. (5): Ratio of integrated flux in fixed velocity intervals eorresponding to the regions of variability in the red and blue 
sides of the rms profile. Col. (6): Full width at quarter- maximum in units of 10^ km s~^. Col. (7): Quarter-maximum shift measured from 
central velocity at quarter maximum in units of 10^ km s~^. Col. (8): Flux-weighted centroid at quarter-maximum in units of 10^ km 
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TABLE 13 
Mkn 668: Broad Ha Properties 



UT Date 


Flux 


V(R) 


P(R) 


F(R)/F(B) 


FWQM 


QMS 


QMC 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1QSQ Till 0*^ 


22itl 


2.7 


n 1 91 -1-0 nofi 


2.16 


11.6 


—0.5 


0.4 




14.2it0.8 


2.7 


n nqO-l-O 00"^ 


2.48 


11.0 


—0.3 


0.6 


iqqn Mav Sf) 


13 7-1-0.7 


2.6 


u . o y _i_ "J . yjvj'j 


2.50 


10.8 


—0.3 


0.6 


IQQp; Ton 91 




2 4 


n 0fi9-l-0 00*^ 


1 ^q 


1 1 n 


— u.o 


n 9 


1 QQf^ ATar 94 


12.8i0.7 


2.0 


n 070-1-0 nOzL 


1.65 


11.0 


—0.1 


0.4 


-Lyyu o Liii 


11.2±0.6 


2.0 




1.53 


11.2 


—0.4 


0.2 


1 QQfi Voh 1 zl 
-LyyvJ j-L^u ±'-t 


14.1it0.8 


1.4 


fl 077-1-0 nOzL 
U.U ( / nzu.uu'i 


1.48 


11.2 


—0.6 


0.0 


1 QQ7 Foh 07 
±yy / p c L) u 1 


7.3it0.4 


1.3 


fl 041 4-0 009 

U . U 'I -L !X; u . u u 


1.48 


11.7 


—0.4 


0.0 


1QQ7 Tun flQ 
xyy ( J nil uy 


8.8it0.5 


1.1 


fl 0^^9-1-0 00*^ 
U.UUZr^U. uuo 


1.48 


10.8 


—0.5 


0.0 


1998 Jan 29 


6.0±0.3 


1.8 


0.034±0.002 


1.46 


11.0 


-0.3 


0.1 


1998 Apr 08 


7.6±0.4 


1.3 


0.040±0.002 


1.40 


11.6 


-0.4 


-0.1 


1998 May 19 


4.7d=0.3 


0.8 


0.024±0.001 


1.46 


11.8 


0.0 


0.2 


1998 Jun 27 


8.6±0.5 


0.5 


0.048±0.003 


1.41 


10.9 


-0.3 


0.1 


1999 Jun 15 


6.0±0.3 


1.2 


0.034±0.002 


1.44 


11.2 


-0.3 


0.1 


2000 Jun 03 


6.3±0.4 


1.8 


0.031±0.002 


1.48 


12.2 


-0.2 


0.1 


2001 Jul 21 


24±1 


1.1 


0.130±0.007 


1.81 


11.4 


-0.3 


0.2 


2002 Jan 08 


29±2 


1.3 


0.154±0.008 


2.17 


11.9 


0.0 


0.5 


2003 Jan 26 


15.9±0.8 


1.6 


0.088±0.005 


2.30 


12.6 


0.3 


0.8 


2004 Mar 29 


21±1 


1.5 


0.111±0.006 


2.25 


12.1 


-0.1 


0.5 


2004 Jun 09 


17.8±0.9 


2.0 


0.093±0.005 


2.34 


12.4 


-0.1 


0.5 



Note. — Col. (1): Date of observation. Col. (2): Flux of broad Ha relative to narrow Ha with la errors. Col. (3): Flux-weighted 
velocity centroid of top 10% of red peak in units of 10^ km s~^. Col. (4): Average flux in top 10% of red peaJc relative to narrow Ha with 
1(7 errors. Col. (5): Ratio of integrated flux in fixed velocity intervals corresponding to regions of strong variability on the red and blue 
side determined from the rms proifle. Col. (6): Pull width at quarter- maximum in units of 10^ km s"!. Col. (7): Quarter -maximum shift 
measured from central velocity at quarter maximum in units of 10^ km s~^. Col. (8): Flux-weighted centroid at quarter-maximum in units 
of 10^ km s-i 



TABLE 14 
3C 227: Broad Ha Properties 



UT Date 


Flux 


V(B) 


P(B) 


F(R)/F(B) 


FWQM 


QMS 


QMC 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1990 Feb 22 


17.8±0.8 


-2.1 


0.120±0.005 


0.25 


13.1 


-0.3 


-0.8 


1993 Dec 13 


22±1 


-2.0 


0.145±0.006 


0.35 


12.3 


-1.2 


-1.1 


1995 Jan 23 


9.8±0.5 


-2.1 


0.067±0.003 


0.33 


11.6 


-1.6 


-1.3 


1995 Mar 24 


14.3±0.6 


-2.1 


0.097±0.004 


0.33 


11.1 


-1.3 


-1.1 


1996 Feb 15 


11.4±0.5 


-2.1 


0.074±0.003 


0.36 


12.4 


-1.5 


-1.3 


1997 Feb 06 


10.6±0.5 


-2.1 


0.072±0.003 


0.36 


12.8 


-1.0 


-1.1 


1998 Jan 29 


15.5±0.7 


-1.8 


0.096±0.004 


0.42 


12.8 


-1.0 


-0.8 


1998 May 20 


7.9±0.4 


-2.0 


0.056±0.003 


0.37 


13.4 


-0.9 


-1.0 


1998 Dec 20 


14.0±0.6 


-1.6 


0.078±0.003 


0.44 


14.4 


-0.7 


-0.6 


1999 Dec 04 


19.1±0.8 


-1.3 


0.115±0.005 


0.51 


14.4 


-0.5 


-0.4 


2001 Oct 25 


16.6±0.7 


-1.1 


0.103±0.005 


0.45 


14.2 


-0.9 


-0.7 


2002 Jan 08 


15.8±0.7 


-1.1 


0.095±0.004 


0.50 


14.9 


-0.3 


-0.3 


2002 Jan 11 


16.3±0.7 


-1.0 


0.098±0.004 


0.50 


14.0 


-0.8 


-0.6 


2002 Apr 12 


11.6±0.5 


-1.1 


0.072±0.003 


0.50 


13.9 


-0.6 


-0.5 


2002 Dec 27 


17.4±0.8 


-0.9 


0.107±0.005 


0.52 


13.2 


-0.8 


-0.5 


2003 Jan 26 


17.2±0.7 


-1.1 


0.104±0.005 


0.52 


13.8 


-1.0 


-0.7 


2003 Apr 05 


15.2±0.7 


-0.9 


0.090±0.004 


0.53 


13.4 


-0.8 


-0.5 


2003 Oct 31 


10.9it0.5 


-0.6 


0.058±0.003 


0.52 


13.4 


-0.7 


-0.4 


2003 Dec 22 


11.5±0.5 


-1.2 


0.075±0.003 


0.43 


12.7 


-0.8 


-0.6 


2004 Mar 28 


11.1±0.5 


-0.9 


0.067±0.003 


0.51 


12.6 


-0.8 


-0.5 



Note. — Col. (1); Date of observation. Col. (2): Flux of broad Ha relative to narrow Ha with la errors. Col. (3): Flux-weighted 
velocity centroid of top 10% of blue peak in units of 10^ km s~^. Col. (4): Average flux in top 10% of blue peak; relative to [O I]A6300 
with la errors. Col. (5): Ratio of integrated flux in fixed velocity intervals corresponding to regions on the red and blue side of the profile 
corresponding to strong variations in the rms spectrum. Col. (6): Full width at quarter-maximum in units of 10^ km s"!. Col. (7): 
Quarter-maximum shift measured from central velocity at quarter maximum in units of 10^ km s~^. Col. (8): Flux- weighted centroid at 
quarter-maximum in units of 10^ km 



Variability of Double-Peaked Emission Lines 



27 
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Velocity (10' km s") 

Fig. 2. Arp 102B: (a) Mean and (b) rnis profile of the broad Ha: line in units of f,y normalized by the total broad-line flux, with dashed 
lines indicating the velocities of the red and blue peaks of the mean profile, measured from the flux-weighted velocity centroid of the top 
10% of the peaks. 
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FlK. 3. — Arp 102B: (a) Total broad-linc flux calibrated relative to the narrow [O I] A6300 line, (b) full width at quarter-maximum, 
(c) shift at quarter-maximum, (d) ratio of the average flux in fixed intervals in the red (+4900 to +6400 km ) and blue (—4500 to 
—5900 km s^^ ) peaks, (e) flux- weighted eentroid at quarter-maximum, (f) fiux- weighted velocity centroid and (g) dispersion of top 10% 
of the blue peak of the broad Hq profiles over the entire duration of our monitoring program. Solid lines show the step-like variation of 
the FWQM that are reproduced by model disk profiles. The least-squares fit sine curves to the oscillations of F(R)/F(B) are plotted with 
solid lines. The times of a rapid flare in total broad-line flux 1989, and a high-state of total broad-line flux in 1989 during which a large 
amplitude flcire occurred, are plotted with open circles. 
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-10 -5 5 10 15 
Velocity (10' km s") 

Fic. 1. 3C 390.3: (a) Mean and (b) rms profile of the broad Hq line in units of normalized by the total broad-line flux, with dashed 
lines indicating the velocities of the red and blue peaks of the mean profile, measured from the flux-weighted velocity centroid of the top 
10% of the peaks. 
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Fig. 5. - 3C 390.3: (a) Total broad-linc flux calibrated relative to the narrow [O I] A6300 line, (b) full width at quarter-maximum, (c) 
shift at quarter-maximum, (d) ratio of the average flux in fixed intervals in the red (-1-3300 to -1-5900 km s"! ) and blue (-1700 to -4500 
km s^^ ) peaks, (e) flux-weighted centroid at quarter-maximum, and (f) peak velocity separation of the broad Ha profiles over the entire 
duration of our monitoring program. 
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Fig. 6. — 3C 390.3: (a, b) Flux-weighted velocity centroid and (c, d) dispersion of the top 10% of the red and blue peaJjs of the 
Ha profiles over the entire duration of our monitoring program. 
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Fig. 7. 3C 390.3: Normalized integrated flux in fixed intervals in the blue wing [F(BW): -4500 to -12,300 kin s^^ ], blue peak [F(B): 
-1700 to -4500 km s"! ], central peak [F(C): -1700 to +15,000 km s"! ], red peak [F(R) :+3300 to +5900 km s"! ], and red wing 
[F(RW): +5900 to +15,000 km ] of the broad Ha profiles over the entire duration of our monitoring program. 
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Fig. 8.— 3C 390.3: -F(B)/F(R) as defined by Veilleux & Zheng (1991) (VZ91), and the velocity of the blue peak for over 30 yr. VZ91 
data measured from the broad H/3 profile are plotted with open circles, and our data measured from the Ha profiles in our monitoring 
program are plotted with filled circles. 



34 



Gezari et al. 




-10 -5 5 10 15 -10 -5 5 10 15 

Velocity Velocity 
(10' km s") (10' km s") 

Fig. 9. — 3C 390.3: A sequence of broad Ha line spectra, without the narrow emission lines removed, which reveal transient sharp 
features in the blue peak. 
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Fig. 10. 3C 332: (a) Mean and (b) mis profile of the broad Ha line in units of /,y normalized by the total broad-line flux, with dashed 
lines indicating the velocities of the red and blue peaks of the mean profile, measured from the flux-weighted velocity centroid of the top 
10% of the peaks. 
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Fig. 11. — 3C 332: (a) Total broad- line flux calibrated relative to the narrow Ha line, (b) full width at half-maximum, (c) shift at 
half-maximum, (d) ratio of the average flux in fixed intervals in the red (-1-7200 to -|-10,000 km s~^ ) and blue(— 5300 to —8000 km s~^ ) 
peaks, (e) flux-weighted centroid at quarter-maximum, and (f) peak velocity separation of the broad Ha profiles over the entire duration 
of our monitoring program. 
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Fig. 12. — 3C 332: Velocity of the (a) red and (b) blue peaJjs measured from the flux-weighted velocity centroid of the top 10% of the 
peaJjs of the broad Ha profiles over the entire duration of our monitoring program. 
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Fk;. 13. — 3C 332: Strong negative eorrelations, measured from the linear correlation eoeffieient R, between the total broad-line flux 
[F(tot)] and tlie pealc velocity separation [l/(R)-y(B)], the full width at half-maximum (FWHM), and the integrated flux in fixed intervals 
in the blue wing (-8000 to -10,800 Itm s'^ ) and the red wing (-1-10,000 to -|-12,800 Itm s^^ ) relative to F(tot). 
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Fk;. 14. — PKS 0235+023: (a) Mean and (b) rms profile of the broad Ha fine in units of /i/ normalized by the total broad-linc flux, with 
dashed lines indicating the velocities of the red and blue peaks of the mean profile, measured from the fiux-weighted velocity centroid of 
the top 10% of the peaJfs. 
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Fig. 15. — PKS 0235+023: (a) Total broad-linc flux calibrated relative to the narrow Hq line, (b) full width at quarter-maximum, (c) 
shift at quarter-maximum, (d) ratio of the integrated flux in fixed intervals in the red (+2300 to +5200 km s"'^ ) and blue(— 3900 to — 
6900 km s^^ ) peaks, (e) fiux-wcightcd centroid at quarter-maximum, (f) fiux-wcightcd velocity centroid and (g) dispersion of top 10% of 
the blue peak of the broad Ha profiles over the entire duration of our monitoring program. 
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FlH. 16. Mkii 668: (a) Mean and (b) rms profile of the broad Hq line in units of /i/ normalized by the total broad-line flux, with a 
dashed line indicating the velocity of the central red peak of the mean profile, measured from the flux-weighted velocity centroid of the top 
10% of the peak. 
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Fig. 17. Mkn 668: (a) Total broad-linc flux calibrated relative to the narrow Ha line, (b) full width at quarter-maximum, (c) shift 
at quarter-maximum, (d) ratio of the integrated flux in fixed intervals in the red side (+3200 to -1-5600 km s^-"^ ) and blue side (—800 to 
—3200 km s^^ ) of the profile, (e) flux-weighted eentroid at quarter-maximum, (f) flux- weighted velocity centroid of top 10% of the central 
red peak of the broad Ha profiles over the entire duration of our monitoring program. 
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Fig. 18. — Mkn 668: Historical velocity curve of the red peak. Data measured for the H/3 profile in 1974 from Osterbrock & Cohen 
(1979), in 1982 from Gaskell (1983), and from 1985 to 1991 from Marziani et al. (1993) are plotted with open circles. Our data measured 
from the Ha profiles in our monitoring program plotted with filled circles. 
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Flc;. 19. — 3C 227: (a) Mean and (b) rms profile of the broad Ha line in units of normalized by the total broad-line flux, with a 
dashed line indicating the velocity of the central blue peak of the mean profile, measured from the flux-weighted velocity centroid of the 
top 10% of the peak. 
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Fig. 20. — 3C 227: (a) Total broad-line flux calibrated relative to the narrow Ha line, (b) full width at quarter-maximum, (c) shift at 
quarter-maximum, (d) ratio of the integrated flux in fixed intervals in the red side (+1700 to -1-3400 km s""*^ ) and blue side (—1200 to 
—4400 km s^^ ) of the profile, (e) flux-weighted centroid at quarter-maximum, (f) flux-weighted velocity centroid of top 10% of the central 
blue peak of the broad Ha profiles over the entire duration of our monitoring program. 
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Fig. 21. — 3C 227: Historical velocity curve of the blue peak of the broad Ha line, with data in 1974 from Osterbrock, Koski, & Phillips 
(1976), and in 1982 from Gaskell (1983) plotted with open circles, and data from our monitoring program plotted with filled circles. 
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Fk;. 22. — 3C 382: (a) Mean and (b) mis profile of the broad Ha line in units of fv normalized by the total broad-line flux in the regions 
of the profile not contaminated by the narrow emission lines. 
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Fig. 23. 3C 382: Integrated flux in the red peak [F(R): +2100 to +3800 km s^^ ], red wing [F(RW); +4400 to +6200, +8000 to +15,000 
km s-i ], and far blue wing [F(BW); -7500 to -11,300 km s"! ] normalized by the total broad-line flux [F(tot)], and the F(RW)/F(BW) 
ratio of the broad Ha profiles over the entire duration of our monitoring program. 
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Fig. 24. — Arp 102B: Circular disk models with changes in the inner and outer radius of the disk (fi and go) that match the observed 
changes in FW QM of the profile, while the other disk parameters remain fixed [i = 31°, g = 3, a = 1050 km s~^). The model parameters 
are defined in ill. 21 
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Fig. 25. — Arp 102B: Velocity curve of the peaks of the model Gaussian excesses oscillating between the red and blue peaks. Solid line 
shows the least-squares fit velocity curve of an orbiting bright spot for the two episodes of oscillation. Both episodes have a best-fit period 
of ~ 2 yr, however the different amplitudes of their velocity curves correspond to a bright spot orbiting at a radius of 427 rgand 757 rg, 
respectively. 
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Fig. 26. — Arp 102B: Integrated flux ratios defined by Sergecv ct al. (2000) measured for tfie broad Ha profiles in our monitoring program. 
Error bars show the mean systematic error in the flux ratios due to the uncertainty in the narrow-line subtractions. Solid lines show the 
bcst-flt sine curves from Sergecv et al. (2000) over the time interval from 1992 to 1995. Dotted lines show the sine curves extended over 
the entire duration of our monitoring program. 
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Fig. 27. — Arp 102B: Change in shape of the profile during the rapid flare in the total broad-line flux in 1989 and a longer flare in 1998. 
Solid line shows the scaled disk profile fit to the broad Ha line. 
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Fig. 28. — Arp 102B: Drift in the blue peak of the broad Ha profile from 1995 to 1996 in comparison to the decrease in peaJc separation 
due to change of outer radius of the disk model from 825 Vgto 1300 rgshown with a thick solid line. The drift of the red peaJc and broadening 
of the peak widths in the disk model are not consistent with the spectra. 
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Fig. 29. — 3C 332: Model circular disk profiles in comparison with the observed broad Ha profiles, with disk parameters of = 190, 
= 540 during the low state, and = 260, = 620 during the high state. 
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Fig. 30. — Mkn 668: Averaged broad Ho profiles that demonstrate the dramatic changes in the red central peak. The thick solid line 
shows the eccentric disk fit to the non-varying portion of the profile. The variations in the profile are modeled with Gaussian excesses 
superposed on the eccentric disk profile, that drift from the red side to the blue side of the profil e and back. These are shown at the bottom 
of the figure, following the line style convention of the legend. See also the discussion in i|4.5.3l 



Variability of Double-Peaked Emission Lines 



55 




Fig. 31. — Mkn 668: Velocity curve of the Gaussian excesses in the eccentric disk profile, in comparison to the velocity curve of a bright 
spot orbiting in the eccentric disk at a radius of 2125 Tg, and with an orbital period of 12 yr. 
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Fig. 32. — 3C 227: Averaged broad Ha profiles that demonstrate the dramatic changes in the blue peak. Gaussian fits to the excess that 
appears to traverse from the blue to the red side of the profile are also shown, following the line style convention of the legend. 
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Fig. 33. — Arp 102B: Broad Ho profiles in units of flux calibrated relative to the narrow [O I] A6300 line, with the subtracted narrow 
lines shown with dotted lines, and the best-fit scaled model disk profile plotted with a thick line. Tick marks show the fiux-weighted 
velocity centroid of the top 10% of the red and blue peaks of the profile. 

APPENDIX 

In this Appendix we present the narrow-hne flux cahbrated spectra of the Ha line profiles of the seven broad-line 
radio galaxies in our spectroscopic monitoring program (Figures 33-52). For all of the objects, except for 3C 382, the 
narrow Ha + [N II]. [O I], and [S II] emission lines have been fitted and removed from the broad Ha lino profiles. The 
subtracted narrow emission lines are shown with a dotted line, and a scaled model for the broad component of the Ha 
line is shown with a solid line. Tick marks indicate the flux-weighted velocity centroids of the peaks of the broad Ha 
profiles. 
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Fig. 34. — Arp 102B: Broad Ha profiles in units of flux calibrated relative to the narrow [O I] A6300 line, with the subtracted narrow 
lines shown with dotted lines, and the best-fit scaled model circular disk profile plotted with a thick line. Tick marks show the fiux-weighted 
velocity centroid of the top 10% of the red and blue peaks of the profile. 
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Fig. 35. — Arp 102B: Broad Ha profiles in units of flux calibrated relative to the narrow [O I] A6300 line, with the subtracted narrow 
lines shown with dotted lines, and the best-fit scaled model disk profile plotted with a thick line. Tick marks show the fiux-weighted 
velocity centroid of the top 10% of the red and blue peaks of the profile. 
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Fig. 36. — Arp 102B: Broad Ha profiles in units of flux calibrated relative to the narrow [O I] A6300 line, with the subtracted narrow 
lines shown with dotted lines, and the best-fit scaled model disk profile plotted with a thick line. Tick marks show the fiux-weighted 
velocity centroid of the top 10% of the red and blue peaks of the profile. 
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Fig. 37. — Arp 102B: Broad Ha profiles in units of flux calibrated relative to the narrow [O I] A6300 line, with the subtracted narrow 
lines shown with dotted lines, and the best-fit scaled model disk profile plotted with a thick line. Tick marks show the fiux-weighted 
velocity centroid of the top 10% of the red and blue peaks of the profile. 
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Fig. 38. — Arp 102B: Broad Ha profiles in units of flux calibrated relative to the narrow [O I] A6300 line, with the subtracted narrow 
lines shown with dotted lines, and the best-fit scaled model disk profile plotted with a thick line. Tick marks show the fiux-weighted 
velocity centroid of the top 10% of the red and blue peaks of the profile. 
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Fig. 39. — Arp 102B: Broad Ha profiles in units of flux calibrated relative to the narrow [O I] A6300 line, with the subtracted narrow 
lines shown with dotted lines, and the best-fit scaled model disk profile plotted with a thick line. Tick marks show the fiux-weighted 
velocity centroid of the top 10% of the red and blue peaks of the profile. 
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Fig. 40. — 3C 390.3: Broad Ha profiles in units of /j^, flux calibrated relative to the narrow [O I] A6300 line, with the subtracted narrow 
lines shown with dotted lines, and the scaled mean profile plotted with a thick line. Tick marks show the flux-weighted velocity centroid 
of the top 10% of the red and blue peaks of the profile. 
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Flc;. 41. — 3C 390.3: Broad Ha profiles in units of /,/, flux caUbratcd relative to the narrow lines, with the subtracted narrow lines shown 
with dotted lines, and the scaled mean profile plotted with a thick line. Tick marks show the flux-weighted velocity centroid of the top 
10% of the red and blue peaJss of the profile. 
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Fig. 42. — 3C 390.3: Broad Ha profiles in units of /j^, flux calibrated relative to the narrow [O I] A6300 line, with the subtracted narrow 
lines shown with dotted lines, and the scaled mean profile plotted with a thick line. Tick marks show the flux-weighted velocity centroid 
of the top 10% of the red and blue peaks of the profile. 
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Fig. 43. — 3C 332: Broad Ha profiles in units of flux calibrated relative to the narrow Ha line, with the subtracted narrow lines 
shown with dotted lines, and the best-fit scaled model circular disk profile plotted with a thick line. Tick marks show the flux- weighted 
velocity centroid of the top 10% of the red and blue peaks of the profile. 
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Fig. 44. — 3C 332: Broad Ha profiles in units of flux calibrated relative to the narrow Ha line, with the subtracted narrow lines 
shown with dotted lines, and the best-fit scaled model circular disk profile plotted with a thick line. Tick marks show the flux- weighted 
velocity centroid of the top 10% of the red and blue peaks of the profile. 
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Fic;. 45. — PKS 0235+023: Broad Ha profiles in units of fi,, flux calibrated relative to the narrow Ha line, with the subtracted narrow 
lines shown with dotted lines, and the best-fit scaled model circular disk profile plotted with a thick line. Tick marks show the fiux-weighted 
velocity centroid of the top 10% of the red and blue peaks of the profile. 
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Fig. 46. — Mkn 668: Broad Hq profiles in units of fu, flux calibrated relative to the narrow Ha line, with the subtracted narrow lines 
shown with dotted lines, and the best-fit scaled model eccentric disk profile plotted with a thick line. Tick marks show the flux-weighted 
velocity centroid of the top 10% of the central red peaJc of the profile. 
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Fig. 47. — Mkn 668: Broad Hq profiles in units of f^, flux calibrated relative to the narrow Ha line, with the subtracted narrow lines 
shown with dotted lines, and the best-fit scaled model eccentric disk profile plotted with a thick line. Tick marks show the flux-weighted 
velocity centroid of the top 10% of the central red peaJc of the profile. 
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Fig. 48. — 3C 227: Broad Ha profiles in units of fi,, flux calibrated relative to the narrow Ha line, with the subtracted narrow lines 
shown with dotted lines, and the scaled mean profile from 1999-2004 plotted plotted with a thick line. Tick marks show the fiux-weighted 
velocity centroid of the top 10% of the central blue peak of the profile. 
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Fig. 49. — 3C 227: Broad Ha profiles in units of fi,, flux calibrated relative to the narrow Ha line, with the subtracted narrow lines 
shown with dotted lines, and the scaled mean profile from 1999-2004 plotted plotted with a thick line. Tick marks show the fiux-weighted 
velocity centroid of the top 10% of the central blue peak of the profile. 
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Fig. 50. — 3C 382: Broad Ha profiles including the narrow lines in units of fi^, normalized by the total broad-line flux in the regions of 
the profile not contaminated by the narrow emission lines, and the best-fit scaled model disk profile plotted with a thick line. 
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Fig. 51. — 3C 382: Broad Ha profiles including the narrow lines, normalized by the broad-line flux, and the best-flt scaled model eccentric 
disk profile plotted with a thick line. 
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Fig. 52. — 3C 382: Broad Ha profiles including the narrow lines, normalized by the broad-line fiux, and the best-fit scaled model eccentric 
disk profile shown in blue. 



